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Electronic  Prototyping  Review  Meeting 

May  23-24,  1995 

Wright  State  University,  Dayton,  Ohio 


May  22: 

6:00  -  8:00  PM  Social  at  Holiday  Inn  in  1st  Floor  Hospitality  Suite 


May  23: 

8:30  Welcome 

Dean  Brandeberry,  Wright  State  University 

8:45  Opening  Remarks:  Electronic  Prototyping 

S.  LeClair.  Wright  Laboratory 

/  \  ■ 

9:00  Soft  Optimization  of  Hard  Problems  (  ^ 

Y.C.  Ho,  Harvard  University 

9:45-  10:00  Break 


10:00  Process  Model  Development  for  Use  with  Discrete  Event  System 
Techniques  ^ 

J.  S.  Gunasekera,  Wright  Laboratory 

10:45  Analyzing  Discrete  Event  Simulation  Models  of  Complex  Manufacturing 
Systems:  A  Computational  Complexity  Approach  ^ 

S.  H.  Jacobson  &  Capt  A.W.  Johnson,  Virginia  Tech 


12:00-1:30  Lunch 


1 :30  Finite  Element  Analysis  for  Simultaneous  (Billet,  Die  and  Part)  Design  for 
Deformation  Processes  y 

T.  Hughes.  Stanford 

2: 15  Automated  3-D  Mesh  Generation  for  Deformation  Process  Design 

M.  Shepherd,  RPI  \/ 

3:00  Alternative  Materials  &  Processes:  Incomplete  Design  and  Multi-process 
Optimization 

A.  Chemaly,  Wright  Laboratory 

3:45-4:00  Break 


4:00  Incomplete  Design  for  Machining  Demonstration 

H.  AL-Kamhawi,  Wright  Laboratory 


May  24: 

8:30  An  Evolving  Computational  Model  for  In-situ  Material  Aging  Assessment 

Y.  Pao,  Wright  Laboratory 

9: 15  Materials  Design  using  Pyramidal  Networks 

A.  Jackson,  Wright  Laboratory 

10:00-10:15  Break 

10:15  A  Rapid  Supervised  Learning  Neural  Network  for  Function 
Interpolation  and  Approximation 

P.  Chen,  Wright  Laboratory 

1 1 :00  Constant-radius  Blending  of  Free-Form  Parametric  Surfaces  and  New 
Results  on  Pythagorean-Hodograph  Curves 

R.  Farouki,  University  of  Michigan 

12:00-1:30  Lunch 

1 :30  A  Constraint-Driven  'Piping- Algorithm'  using  Cyclides 

D.  Dutta,  University  of  Michigan 


2:30-4:00  Forum  Discussion: 

User  Perspectives  on  Electronic  Prototyping  Needs. 

Youngstown  State  University ,  Technology  Development  Consortia  (TDC) 
San-Antonio  Air  Logistic  Center  (SA-ALC) 

Sacramento  Air  Logistic  Center  (SM-ALC) 

Oklahoma  City  Air  Logistic  Center  (SA-ALC) 

Systems  Research  Laboratories  (SRL) 

Allison  Engines 

BF  Goodrich  (Wheel  &  Brake  Division) 

Lockheed-Martin 

Cherry  Point  Naval  Depot 

General  Motors  Corporation  -  Saginaw  Division 

Tennessee  Technological  University  -  Center  for  Mfg  Research 

Southern  Illinois  University 


& 


c 

0) 

0) 

0) 


o> 

O) 


c: 

o 

•  ■Ml 

4m«| 

(0 


o 

(G 


0 

(S 

o 

•g 

o 


(0 


(G 

E 

w 

(/) 

jG 

o 

Tg 

o 


S 

cc 

(/) 

O 

LL 

< 


C  +3 
O  <1> 

«  o 

Q_  Ui 

I  o 
8.2 
M-  (0 

o  >* 

(G  (G 
<D  C 

®  c 

:£  <G 

c  O) 
—  c 

o  <i> 
!=  “O 
<0  o 

®  c 
(0  c 

O 

a:  o 


jO 

O 

0) 

u 

o 

> 

0) 

CO 


(G 

O) 

•  Mi 

(0 

0) 


(G 

o| 

c 
■  ■■■ 
u.' 

0. 

Q> 

(G 

V. 

oi 

4i4 

o 

0 

^1 

Q 

(0 

(G 

flHn 

0) 

4— 

<G 


(n 

<D 

^  0  CO 
^  (0 

.E  >o 
o 

0.0*^ 

S  s  “ 

o  Q.(G 

OC  «  *=: 

-j  ® 

0  c 
(0  .5  ? 

0  CO  C 

CD  0  O 

o*  0  E 

0  w  o 

08*5 

0  £  i  I 

^0-  ‘j= 

x:  tf>  00 

0  ■:5  #- 
V.  0  c 

0  'z  n 
0  0  O 
0  «  ‘*- 
0  0  ««- 

QC  3  o 

0  c 
=  0-0 
Q.0  g; 
^0  0 


■O 

0 

■O 


0 

»_ 

(D 

CC 

CO 

O 

u. 

<1 

c 

0 

0 

n 

0 

> 

0 


0 

0 
•  ■MB 

4i4 

0 

0 

>: 


O 

0 

H 

0 

c 

*5) 

mmm 

> 

eS 


0 

> 

0 


JZ 

o 


0 

3 

■O 

3 

CL 


0 

C 

o 

o 


University  Visiting  Scientist  Participants  who  have  been  WL  FuncJed: 
Case  Western  Reserve,  Wright  State,  Ohio  State,  UCinn. 
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design,  manufacture  and  repair  our  defense  systems, 

we  must  also  find  ways  to  automate  the  organization  of 

new  experience(s). 
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Optimized  tn-situ  inspection 
Probe  Path  Pianning 
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Discrete  Optimization 
Method 


State-Space”  Process  Design  And  Control  Optimization 

Methodology  Verification 


Self-Directed  Control/Process  Discovery  for  Thin-film  Processing  -  using  FLN 
technology  to  automate  non-linear  regression  and  address:  1 )  error  characterization, 

2)  in-situ  adaptation,  and  3)  problem  chai  aclerization/discovery.  The  essence  of  these 
contributions  will  be  to  automate  materials  research  and  augment  researchers  in 
improving  in-house  research  productivity  and  affordability  of  new  Materials  &  Processes. 


The  First  to  the  Next  75  Years 


Air  Force  Wright  Laboratory 
Materials  Directorate 
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Controlled  Composition  &  Morphology 

Integrated  (2D&  3D)  Nano  -  Macro  Scale 


Air  Force  Patent  Inventions 
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Soft  Optimization  of  Hard  Problems 

by  Y.C.  Ho 
Harvard  University 

Abstract 

It  can  be  argued  that  optimization  in  the  general  sense  of  making  things 
BETTER  is  the  principal  driver  behind  all  of  prescriptive  scientific  and 
engineering  endeavor,  be  it  operations  research,  control  theory,  and 
engineering  design.  It  is  also  true  that  the  real  world  is  full  of  complex 
decision  and  optimization  problems  that  we  cannot  solve.  While  the 
literature  on  optimization  and  decision  making  is  huge,  much  of  the 
concrete  analytical  results  are  associated  with  what  may  be  called  Real- 
Variable-Based  methods.  The  idea  of  successive  approximation  to  an 
optimum  (say,  minimum)  by  sequential  improvements  based  on  local 
information  is  often  captured  by  the  metaphor  of  "skiing  downhill  in  a  fog". 
The  concepts  of  gradient  (slope),  curvature  (valley),  and  trajectories  of 
steepest  descent  (fall  line)  all  require  the  notion  of  derivatives  and  are  based 
on  the  existence  of  a  more  or  less  smooth  response  surface.  There  exist 
various  first  and  second  order  algorithms  of  feasible  directions  for  the 
iterative  determination  of  the  optimum  (minimum)  of  an  arbitrary  multi¬ 
dimensional  response  or  performance  surface.  Considerable  number  of 
major  success  stories  exist  in  this  genre  including  the  Nobel  prize  winning 
work  on  linear  programming. 

On  the  other  hand,  we  submit  that  the  reason  many  real  world 
optimization  problems  remain  unsolved  is  pardy  due  to  the  changing  nature 
of  the  problem  domain  in  recent  years  which  makes  calculus  or  Real- 
Variable-Based  methods  less  applicable,  e.g,, ,  a  large  number  of  human- 
made  system  problems,  such  as  manufacturing  automation,  communication 
networks,  computer  performances,  and/or  general  resource  allocation 
problems,  involve  combinatorics  rather  than  real  analysis,  symbols  rather 
than  variables,  discrete  instead  of  continuous  choices,  and  synthesizing  a 
configuration  rather  than  proportioning  the  design  parameters.  Such 
problems  are  often  characterized  by  the  lack  of  structure,  presence  of  large 
uncertainties,  and  enormously  large  search  space.  Optimization  for  such 
problems  seem  to  call  for  a  general  search  of  the  performance  terrain  or 
response  surface  as  opposed  to  the  "skiing  downhill  in  a  fog"  metaphor  of 
Real-Variable-Based  performance  optimization  (footnote  #1)  .  Arguments 


IS 


can  also  be  made  on  tJie  technological  front.  Sequential  algorithms  were 
often  dictated  as  a  result  of  the  limited  memory  and  centralized  control  of 
earlier  generations  of  computers.  With  the  advent  of  modern  general 
purpose  parallel  computing  and  essentially  unhmited  size  of  virtual 
memory,  distributed  and  parallel  procedures  or  a  network  of  machines  can 
work  hand-in-glove  with  Search-Based  methods  of  performance  evaluation. 

It  is  one  of  the  thesis  of  this  talk  to  argue  for  such  a  complementary 
approach  to  optimization. 

If  we  accept  the  need  for  Search-Based  methods  as  a  complement  to  the 
more  established  Real-Variable-Based  analytical  techniques,  then  we  can 
next  argue  that  to  quickly  narrow  the  search  for  optimum  performance  to  a 
"good  enough"  subset  in  the  design  universe  is  more  important  than  to 
estimate  accurately  the  values  of  the  system  performance  during  the  initial 
stages  of  the  process  of  optimization.  We  should  compare  order  first  and 
estimate  value  second,  i.e.,  ordinal  optimization  comes  before  cardinal 
optimization.  Colloquial  expressions,  such  as  "ballpark  estimate",  "80/20 
solution"  and  "forest  vs.  trees",  state  the  same  sentiment.  Furthermore,  we 
shall  argue  that  our  preoccupation  with  the  "best"  may  be  only  an  ideal  that 
will  always  be  unattainable  or  not  cost  effective.  Real  world  solutions  to  real 
world  problems  will  involve  compromise  for  the  "good  enough"  . 

In  summary,  the  purpose  of  this  talk  is  to  establish  the  distinct  advantages 
of  the  softer  approach  of  ordinal  optimization  for  the  Search-Based  type  of 
problems.  More  specifically,  the  presentation  will  address  general 
properties,  show  the  many  orders  of  magnitude  improvement  in 
computational  efficiency  that  is  possible  under  this  mind  set,  and  finally 
advocate  the  use  of  Soft  Optimization  as  a  complement  to  the  more  traditional 
cardinal  approach  to  engineering  design.  Example  applications  such  as  rare 
event  probability  in  systems,  ATM  switches  in  communication,  buffer 
allocation  in  manufacturing,  and  control  system  design  will  be  highlighted 
to  substantiate  the  case  (footnote  #2).  Relationships  to  other  subjects  such  as 
fuzzy  sets,  genetic  algorithms,  and  rank  statistics  will  also  be  mentioned. 

Footnotes 

#1  We  hasten  to  add  that  we  fully  realize  the  distinction  we  make  here  is  not  absolutely 
black  and  white.  .A  continuum  of  problem  types  exist.  Similarly,  there  is  a  spectrum  of 
the  nature  of  optimization  variables  or  search  space  ranging  from  continuous  to  integer 
to  discrete  to  combinatorial  to  symbolic. 

#2  Joint  work  with  Ph.D.  students,  Mei  Deng,  M.  Larson,  E.T.W.  Lau,  L.H.  Lee,  D.W.G.  Li, 
N.  Patsis 
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Performance  Evaluation 
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Universal  Alignment  Probability  Curves 
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Sampling  &  Search  in  Strategy  Space 
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Process  Model  Development  for  use  with  Discrete 
Event  System  Techniques 

by  Jay  S.  Gunasekera, 

Ohio  University  and 
Visiting  Research  Scientist 
Wright  Laboratory 

Abstract 

In  the  past  decade,  process  analysis  and  modeling  has  enjoyed  a 
resurgence  owing  largely  to  the  continued  growth  in  the  performance-to- 
cost  ratio  of  computational  systems.  Traditionally,  process  modeling  has 
been  limited  to  analysis  of  a  single  step  in  a  typical  material  processing 
sequence,  and  process  optimization  has  generally  been  carried  out 
empirically  by  iteratively  testing  and  improving  designs  for  each  step  in 
the  process. 

The  current  investigation  is  aimed  at  developing  process  designs  which 
will  enable  the  analysis  and  optimization  of  a  material  processing  sequence. 
Typically,  such  a  sequence  may  involve  several  stages  and  alternative 
processing  paths  for  a  given  component.  It  may  also  involve  both  Real- 
Variable  and  Search-Based  (discrete  event)  models.  Of  principle  research 
interest  is  that  the  material  processing  sequence  be  globally  optimized  to 
achieve  improvements  in  affordability,  producibility  and  performance  of 
Air  Force  Systems. 

Material  process  models  are  capable  of  representing  and  analyzing 
complex  material  flow  problems  (2D  and  3D)  and  non-isothermal  ' 
deformation  processes  with  a  high  degree  of  accuracy.  Material  process 
models  are  generally  non-linear  and  therefore,  rely  on  numerical 
techniques  for  solutions.  These  models  attempt  to  simulate  the  kinetics  of 
material  processing,  and  typically  consist  of  inputs,  outputs,  constraints  and 
a  set  of  equations  to  be  differentiated  and/or  integrated  which  describe  the 
physics  of  the  process.  The  model  may  be  treated  as  a  "Black  Box"  by  the 
optimization  algorithm.  The  entire  manufacturing  process  for  a  "state-of- 
the-art"  engine  part  (ex:  Integrated  Blade  Rotor,  IBR)  consists  of  several 
individual  discrete  processes,  such  as  casting,  HIPing,  forging,  heat  treating, 
machining  etc.  Each  of  these  processes  has  specific  inputs,  outputs, 
equipment  settings,  constraints  and  workpiece  information.  The  model 


incorporates  all  of  this  information  internally  and  generates  variables  for  a 
set  of  specific  inputs. 

There  are  several  methods  available  for  the  analysis  and  modeling  of 

* 

material  deformation  processes.  They  include  lumped  models  (i.e.,  no 
distribution  of  variables),  slab  models  (i.e.,  geometric  simplifications),  upper 
bound  models  (UBET)  models,  finite  difference  models  and  finally,  finite 
element  models.  The  former  models  are  very  fast  in  terms  of  computer 
execution  times  but  provide  less  detailed  information,  whereas  the 
latter  models  are  very  accurate  and  are  capable  of  predicting  the 
variables  in  2D  or  3D,  but  are  very  time  consuming.  The  fundamental 
approach  in  the  discrete  event  simulation  model  is  to  use  models  which  are 
fast  during  the  early  phase  of  optimization  and  progressively  increase  the 
accuracy  of  the  models  when  the  optimum  is  reached. 

This  paper  will  demonstrate  the  development  of  multiple  levels  of  models 
for  forging  and  heat  treatment  processes.  The  models  are  being 
incorporated  into  an  overall  discrete  event  simulation  model  and  will  be  used 
to  determine  the  global  optimum  for  the  entire  material  processing 
sequence  which  will  ultimately  involve  billet  preparation  and  finish  part 
machining  operations. 


PROCESS  MODEL  DEVELOPMENT 
FOR  USE  WITH  DISCRETE  EVENT  TECHNIQUES 


Jay  S.  Gunasekera 
Wright  Laboratory 
Wright-Patterson  AFB 
May  23, 1995 


•  The  overall  goal  of  this  research  is  to  establish  a  systematic  methodology  for 
optimizing  complex  manufacturing  processes  with  respect  to  affordability  and 
quality. 


•  This  research  focuses  on  the  development  of  analytical  and  numerical 
models  for  a  wide  range  of  material  processes  which  are  important  to 
production  of  Air  Force  systems. 


Photo  of  an  extruded-blade  EBR  during  fabrication. 


•  Components  for  Air  Force  systemstypically  involve  advanced  materials  and 
intricate  geometries  with  specified  microstructures  and  mechanical  properties. 

•  Photo  of  integral  blade  and  rotor  component. 

•  Drivers  for  System  Performance 

-  Ability  to  effectively  produce  integral  structure  components 

-  Ability  to  effectively  make  shapes  from  advanced  materials 

-  Ability  to  achieve  Microstructure  &  Property  control  in  “Large 
Sections” 


COMPARISON  OF  OPTIMIZATION  APPROACHES 


New 

Approach 


Traditional 

Approach 


Casting 


Forging 


Heat  Treatment 


Machining 


Global  Optimum 


Traditional  Approach: 

•Start  with  a  feasible  process  as  baseline. 

•Optimize  each  step  in  the  process  individually  with  knowledge  of  the 
processing  operation. 

•Perform  some  trade-offs  between  the  processing  steps  if  of  obvious 
advantage. 

•Hope/Assume/Pray  the  result  is  a  global  optimum  for  the  process,  r 
New  Approach: 

•Treat  each  step  as  a  single  event  in  a  series  required  to  produce  the 
part. 

•Perform  optimization  of  the  entire  series  by  explicitly  trading  off  the 
individual  steps. 


NOMINAL  IBR  PROCESSING  STEPS 


Typical  EBR  forming  operations  with  associated  shapes. 

•Simple  IBR  shape  from  many  processing  steps. 

•Each  step  can  be  broken  down  into  events. 

•Each  event  can  be  further  broken  down  into  multiple  material  behavior  events 
that  determine  flow  (pattern  filling)  and  microstructure. 

•On  this  level. ..the  entire  process  is  a  series  of  events  that  occur  to  produce  the 
final  part. 

•The  goal  is  to  optimize  the  process. 


ACCURRACY 


Accuracy  Vs. 
Execution  Time  of 
Various  Models 
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I  Updated 
I  UBET 


I  Slab 


Distributed  | 
Lumped  | 


I  Lumped 


EXECUTION  TIME 


Modeling  Methods  for 
Deformation 
Processes _ 


•  Slab 

-  No  Motion 

-  Predicts  Die  Forces 

-  Cannot  Predict  Fill 

•  Updated  UBET 

-  2  or  3  Dimensional  Models 

-  Capable  of  Rough  Property 
Distribution  Predictions 

-  Capable  of  Simulating  Die 
Movement  and  Predicting 
Geometry  and  Fill 

•  FEM 


PLASTICITY 
Governing  Equations 

Equilibrium 

da,. 

—^=0 
dx, 

Yield  Criterion 

Constitutive  Equations 


Compatibility  Conditions 


MODIFIED  SLAB  ANALYSIS 
Governing  Equations 


Vertical  Equilibrium 
-Cj  =  P;  (l-tan  a) 

=  P„  (l+tan  P) 
Horizontal  Equilibrium 


+  -a^  +  +  P„  tan  p  -  t„-P 

ar  r  ar 


—  =  (tan  p-tan  a) 
dr 


Coulomb  Friction 


Sticking  Friction 
=  mk 


^  tan  a-  t^=0 


UPDATED  UPPER  BOUND  ELEMENTAL 

TECHNIQUE  (UBET) 

UBET  Energy  Minimization 

UBET  Velocity  Field  Elements 


UBET  Velocity  Fields 
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UBET  Energy  Minimization 
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w.  =  Internal  Plastic  Deformation  Power 
w^  =  Shear  Losses  Along  Element  Boundaries 
viy  =  Friction  Losses  at  Die -Workpiece  Interface 


Compatibility 


Program  Structure 


LINE(P14,P15) 


Targeted  Schedule 
for  Process  Model  Development 
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Analyzing  Discrete  Event  Simulation  Models  of 
Complex  Manufacturing  Systems: 
a  Computational  Complexity  Approach 

by  Sheldon  H.  Jacobson  and  Alan  W.  Johnson 
Virginia  Tech 

Abstract 

Discrete-event  computer  simulation  is  a  powerful  modeling  tool  for 
analyzing  complex  systems  that  cannot  be  studied  analytically.  The  value  of 
simulation  depends  on  how  effectively  and  efficiently  one  can  build  and 
utilize  the  simulation  model  to  gain  insights  into  complex,  real-world 
systems. 

The  goal  of  this  research  is  to  introduce  a  new  set  of  tools  for  building  and 
analyzing  simulation  models  of  complex  systems,  with  a  particular  emphasis 
on  manufacturing  processes  for  turbine  engine  parts,  under  development  at 
Wrighf  Laboratory  (Materials  Process  Design  Branch).  The  basic  research 
question  to  be  considered  is:  Given  a  corriplex  manufacturing  system,  where 
the  different  processes  and  operations  can  be  modeled  using  discrete  event 
simulation  model  events,  is  it  possible  to  determine  the  sequence  of  processes 
and  operations,  hence  events,  that  results  in  finished  products  with 
prespecified  attributes  and  output  measures?  The  attributes  and  measures  of 
such  manufacturing  processes  typically  involve  physical  and  operational 
properties  of  the  finished  products,  as  well  as  metrics  for  affordability, 
predictability,  and  repeatability  Of  the  outputs. 

Research  questions  of  this  type  are  often  addressed  using  adhoc 
procedures  and  problem-specific  approaches.  This  project  uses  the  theory  of 
computational  comple.xity  to  provide  a  framework  to  study  these  research 
questions,  as  well  as  construct  new  heuristic  algorithms  to  address  them. 
These  heuristics  involve  simulated  annealing,  genetic  algorithms,  and  tabu 
search.  These  three  heuristics  will  be  cross-fertilized,  resulting  in  hybrid 
approaches,  as  well  as  enhanced  to  exploit  any  special  structure  offered  by 
the  manufacturing  problems. 
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BACKGROUND  AND  MOTIVATION 


*  RT  ITT  nTNG  AND  ANALYZING  DISCRETE  EVENT 
SIMULATION  MODELS  ARE  DIFHCULT  TASKS. 


*  COMPUTER  SUPPORT  (SOFTWARE)  AND 
COMPUTING  POWER  (HARDWARE)  ARE 
PLENTIFUL. 


*  RESEARCHERS  HAVE  BEEN  WORKING  ON 

DEVELOPING  A  COMPREHENSIVE  MODELING 
FRAMEWORK  TO  AID  PRACTmONERS  BUILD 
AND  ANALYZE  SIMULATION  MODELS. 
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DEFINITIONS 


*  A  SYSTEM  IS  A  COLLECTION  OF  INTERACTING 

ENTITIES  DETERMINED  BY  CERTAIN  LAWS 
AND  POLICIES. 

*  THE  STATE  IS  A  COMPLETE  DESCRUTTON  OF  A 

SYSTEM. 

*  AN  EVENT  IS  WHERE: 

-  STATE  CHANGES  OCCUR 

-  OTHER  EVENTS  ARE  SCHEDULED/CANCELLED. 


*  A  MODEL  SPECIFICATION  (MS)  IS  A  SYSTEM 
REPRESENTATION  OF  A  DISCRETE  EVENT 
SIMULATION  MODEL. 


*  A  MODEL  IMPLEMENTATION  (MI)  IS  A  COMPUTER 
REPRESENTATION  OF  A  MODEL 
SPECIFICATION. 


ArrRSSTBILITY: 


Tn  stance  I  A  discrete  event  simulation  MS  with 

associated  MI 
An  initial  event  Eo 
A  state  S 

A  non-negative  finite  integer  M. 

Problem:  Find  a  sequence  of  events  Ei,E2v..,Em,  m^ 

such  that  the  execution  of  the  sequence 
yields  E0E1E2  •  •  •  Em  S. 


*  TfflS  IS  A  TRAJECTORY  PROBLEM. 

*  WHY  IS  THIS  PROBLEM  INTERESTING? 


MANUFACTURING  APPLICATION 
TURBINE  engine  MANUFACTURING 

*  BILLETS  ARE  TRANSFORMED  INTO  FINISHED 

PRODUCTS. 

*  TITANIUM  ALLOYS  (TiAI)  ARE  USED  FOR  THESE 

PRODUCTS. 

*  SEVERAL  DISCRETE  STEPS  (EVENTS)  ARE 

INVOLVED  IN  THIS  TRANSFORMATION: 

-  PRECONDITIONING 

-  TEMPERATURE  VARIATIONS 

(HEAT  TREATMENTS,  COOLING) 

-  FORGING  AND  EXTRUSION 


-  997 


FUNDAMENTAL  RESEARCH  QUESTION: 

*  HOW  CAN  THESE  STEPS  (EVENTS)  BE  SEQUENCED 
SUCH  THAT  THE  FINISHED  PRODUCT  MEETS 
CERTAIN  GEOMETRIC,  AND  MICRO- 
STRUCTURAL  PROPERTIES? 

-  AFFORDABILITY 

-  PREDICTABILITY 

-  REPRODUCIBILrTY 


SOLUTION: 

*  FORMULATE  THE  PROBLEM  AS  A  (DYNAMIC) 

COMBINATORIAL  OPTIMIZATION  PROBLEM, 
WHERE  THE  DISCRETE  OBJECTS  ARE 
SIMULATED  EVENTS. 

*  DEVELOP  ALGORITHMS  THAT  WILL  ADDRESS  THIS 

(SEQUENCING)  PROBLEM  (ACCESSIBILITY). 


Steepest  Descent  Algorithm 

(Mmimization  Problem) 


-  Select  an  initial  incumbent  solution  a  from  solution  space  E  (o-  €  1) 

-  Define  a  cost  function  f((r)  for  solution  a 

-  Select  a  solution  neighborhood  ftmction  that  defines  all 

neighbors  a*  of  a  (a*  € 

-  Select  a  stopping  criterion,  N,  that  defines  the  total  number 

of  iterations  executed 


Set  the  repetition  counter  n  =  0 
While  n  7^  N: 

Randomly  generate  a  solution  (f  G  rj((j) 
Calculate  cost  function  change  8  =  f{(f)  -  f((T) 
if  8  <  0,  then 
else  n«=n+l 


Source:  Eglese,  R.W.  "Simulated  Annealing:  A  Tool  for  Operationai  Research",  European  Jotmtal  of 
Operational  Research,  Vol  46,  pp  271-281,  1990. 


Tabu  Search  Algorilhrii 

(Minimization  Problem) 


-  Select  an  initial  incumbent  solution  <r  from  solution  space  E  (<r  €  E) 

-  Define  optimal  solution  as  a* 

-  Define  a  cost  function  f(a)  for  solution  (r 

-  Select  a  solution  neighborhood  function  that  defines  all 

neighbors  a*  of  (t  (<r^  €  jj((r)) 

-  Define  a  list  L((r)  of  length  k  whose  elements  are  called  Uil>u  moves 


1 .  -  Set  the  repetition  counter  n  =  0 

-  Initialize  L(c7)  as  empty 

-  <7*  ^  O’ 

2.  - //t7(o)  -  L((7)  =  0  go  to  step  4. 

-  else 

—  n  <=  n+1 

-  Select  (fa  €  77(0-)  -  L{(t)  such  that 

f((fn)  =  minimum {f((7*):  (f  €  -  L((y)} 

3. 

-  if  f(o)  <  then  &^^<r 

4.  -  z/stopping  criterion  met,  or 

if  77(0)  -  L((7)  =  0  (from  step  2.),  then  stop 

-  else  update  tabu  list  L((t),  and  return  to  step  2. 

Source:  Glover,  F.,  "Tabu  Search  -  Part  1",  ORSA  Journal  on  Computing,  Vol  1,  No.  3,  pp  190-206,  1989 


KEY  COMPONENTS  OF  ALGORITHMS 


*  CONFIGURAHON  SPACE  APPROACH 


-  SOLUTION  SPACE 

-  SET  OF  ALL  POSSIBLE  SEQUENCES  OF  EVENTS 
(LARGE!) 


-  COST  (OBJECTIVE)  FUNCTIONS 

-  MUST  CAPTURE:  -  ALL  GEOMETRIC  AND 

MICROSTRUCTURAL 
PROPERTIES  OF  THE 
FINISHED  PRODUCT 
-  KEY  CHARACTERISTICS 
OF  THE  MANU¬ 
FACTURING  PROCESS 


-  NEIGHBORING  SOLUTIONS 
-  IDENTIFY  MECHANISMS  TO  EITICIENTLY 

traverse  the  solution  space 


*  algorithm  research  issues 

-  VARIANCE  ENHANCEMENT  TECHNIQUES 
&  PROBLEM  RANDOMIZATION 

-TUNNELLING 

(RANDOM  PERTURBATION  RELAXATION) 

-  SENSITIVITY  ANALYSIS  OF  ALGORITHMS 

-  FALSE  NEGATIVE  PROBLEM 
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Finite  Element  Analysis  for  Simultaneous 
(Billet,  Die  and  Part) 

Design,  for  Deformation  Processes 

Thomas  J.R.  Hughes 
Division  of  Applied  Mechanics 
Durand  Building 
Stanford  University 
Stanford,  California  94305 

Abstract 

This  presentation  describes  a  research  program  directed  towards  the 
development  of  automated  design  procedures  for  aluminum  extrusion 
technology.  The  objective  is  to  eliminate  costly  trial  and  error  by  being  able  to 
simultaneously  design  the  product,  die,  billet,  and  process  (e.g.,  extrusion 
temperatures  and  speeds,  uniformizing  metal  flow,  etc.),  within  constraints  of 
feasibility,  and  satisfying  objectives  including,  but  not  limited  to,  optimizing 
shape,  surface  finish,  and  properties  of  the  product,  processing  costs,  time  to 
market,  and  full  utilization  of  capabilities.  The  approach  is  based  on  the 
development  of  efficient  and  effective  analysis  of  the  whole  processing  system 
employing  newly  developed  finite  element  solution  technologies  for  complex, 
multiregion,  multiphysical  behavior.  Generalizations  of  these  methodologies  to 
include  Arbitrary  Lagrangian-Eulerian  (ALE)  mesh  descriptions  for  nonlinear, 
elasto-viscoplastic  mechanical  constitution  equations  will  allow  the  faithful 
modeling  of  the  metal  flow  within  the  die  system  and  the  accurate  attainment  of 
final  shape  upon  exit.  Automatic  meshing  and  adaptive  remeshing  will  insure 
efficient  and  accurate  simulation  of  the  entire  forming  process.  New  element 
technologies  facilitating  the  use  of  general  meshing  procedures  for  difficult 
metal-forming  processes  involving  a  variety  of  kinematical  constraints,  such 
as  incompressibility,  contact,  etc.,  are  utilized.  Feature  based  design 
methodologies,  parametric  modeling,  and  knowledge-based  engineering 
techniques  will  constitute  the  fundamental  methodologies  for  representing 
designs,  managing  the  hierarchy  of  analysis  models,  performing  model 
reduction  and  feature  removal,  and  effectively  utilizing  design  knowledge. 
Modern,  three-dimensional,  interactive  visualization  procedures  are  employed 
to  animate  simulations  and  design  evolution.  The  software  will  be  accessed 
through  a  modern  and  easy-to-use  graphical  user  interface  developed  for 
extruders,  die  makers,  and  designers.  The  current  status  of  our  software  and 
plans  for  executing  the  research  and  development  program  will  be  presented. 


Finite  Element  Analysis  for  Simultaneous 
(Billet,  Die  and  Part)  Design 
for  Deformation  Processes 


T.  Hughes 
Stanford  University 


Outline 


•  Introduction 


•  Technical  Discussion 


•  Singularities 

•  Localization  of  Deformation 


•  Adaptive  Refinement 

•  Element  Technology 

•  Contact 


•  Friction 


•  ALE 


Outline  (cont.) 


•  Technical  Discussion  (cont.) 

•  Space-time  Formulation 

•  Die  Flexibility 

•  Region-based  Computational  Architecture 

•  Parallel  Processing 

•  Equation  Solving  Technologies 

•  Visualization 


•  Control  of  Process 


Outline  (cont.) 


•  Leveraging  Existing  Technology 
•  Modeling  Environment 


•  Mesh  Generation  and  Adaptivity 

•  Geometric  Modeling 

•  Simulation  and  Visualization 


•  User  Interface 


•  Optimization 


•  Summary 


Introduction 


•  Research  directed  towards  the  development  of  automated 
design  procedures  for  aluminum  extrusion  technology 

•  Goal:  Eliminate  costly  trial  and  error  by  being  able  to 
simultaneously  design  the  product,  die,  billet,  and  process, 
within  constraints  of  feasibility  and  satisfying  objectives 
including  but  not  limited  to  optimizing  shape,  surface 
finish,  and  other  product  properties,  processing  costs,  time 
to  market,  full  utilization  of  capabilities,  etc. 

•  Approach  is  based  on  the  efficient  and  effective  analysis 
of  the  whole  process  employing  newly  developed  finite 
element  solution  technologies  for  complex,  multiregion, 
multiphysical  behavior 


Introduction  (cont.) 


•  ALE  mesh  descriptions  for  nonlinear,  elasto- viscoplastic 
mechanical  constitutive  equations  will  allow  accurate  mod¬ 
eling  of  the  metal  flow  within  the  die  and  container  and  the 
accurate  attainment  of  final  shape  upon  exit 

•  New  element  technologies  facilitating  the  use  of  general 
meshing  procedures  for  difficult  metal-forming  processes 
involving  kinematical  constraints,  such  as,  incompressibility 
and  contact,  are  utilized 

•  Automatic  meshing  and  adaptive  remeshing  will  be 
employed  to  attain  efficient  and  accurate  simulations  (see 
presentation  of  M.  Shephard) 


Introduction  (cont.) 


•  Feature  based  design  methodologies,  parametric  modeling, 
and  knowledge  based  engineering  technologies  will  be  uti¬ 
lized  to  represent  designs,  manage  hierarchies  of  analysis 
models,  perform  model  reduction  and  feature  removal,  and 
effectively  utilize  experience  and  design  knowledge  (see 
presentation  of  A.  Chemaly) 

•  Modern,  three-dimensional,  interactive  visualization 
procedures  will  be  employed  to  animate  simulations  and 
design  evolution 

•  The  software  will  be  accessed  through  an  easy  to  use 
graphical  user  interface  developed  for  extruders,  die  makers 
and  designers 


Technical  Discussion 


Singularities 


•  Geometrical  discontinuities  within  die  and  at  the  point 
of  exit  give  rise  to  “attached”  singularities  in  the  rate-of- 
deformation  field 

•  These  have  the  structure  of  classical  elliptic  singularities 

•  Approach:  p- adaptivity 
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Localization  of  Deformation 


•  In  fiat-faced  dies,  very  strong,  localized  gradients  exist 
across  “shear  surfaces,”  causing  a  “dead-metal”  zone  to 
form  in  front  of  the  die 

•  Classical  polynomial  interpolation  is  not  well-suited  to 
approximate  such  abrupt  discontinuities 

•  Special  ingredients  need  to  be  incorporated  in  a  numerical 
formulation  to  avoid  mesh-dependent  pathologies 

•  Various  methodologies  have  been  proposed 

•  Approach:  Combination  of  /i-adaptivity  and  subgrid- 
scale  modeling  techniques  which  embody  the  features  nec¬ 
essary  to  avoid  mesh  sensitivity 


/ip- Adaptive  Refinement 


•  Both  h  and  p  strategies  are  relevant  in  simulating  extrusion 
processes 

•  /i- adaptivity  for  interior  surfaces  of  near-discontinuity 

•  p-adaptivity  for  elliptic  singularities  and  high-accuracy  in 
regions  of  smooth  response 

•  Approach:  Combined  /ip-adaptive  strategy 


Element  Technology 


•  Incompressible  kinematics  is  fundamental  to  the  description 
of  the  billet  and  part 

•  Fully  automatic  meshing  and  adaptive  remeshing  are  avail¬ 
able  for  tetrahedra 

•  Need  for  /ip-adaptive  tetrahedral  meshing  puts  stringent 
requirements  on  element  technology  for  extrusion  problems 

•  Most  combinations  of  velocity  and  pressure  interpolations 
give  rise  to  unstable,  hence  nonconvergent,  numerical 
approximations  within  the  classical  Galerkin  formulation 

•  Approach:  Stabilized  methods  which  support  any  com¬ 
bination  of  velocity  and  (continuous  or  discontinuous) 
pressure  (H.-Franca-Stenberg-Brezzi-etc.) 


Contact 


•  Mechanical  and  thermal  contact  between  billet  and 
container  and  die  plays  a  significant  role  in  the  extrusion 
process 

•  The  mathematical  sensitivities  are  similar  to  that  for 
incompressible  behavior 

•  The  standard  treatment  of  contact  within  the  Galerkin 
method  does  not  provide  a  suitable  basis  for  /ip-strategies 

•  Approach:  Stabilized  methods  (H.-Barbosa) 


Friction 


•  Friction  is  an  important  aspect  of  billet-container-die 
interaction 

•  Friction  is  an  important  source  of  heat  generation 

•  Coulomb  friction  is  viewed  as  inadequate  for  representing 
extrusion  processes 

•  Approach:  Fundamental  study  in  comparing  mechanical 
and  thermal  implications  of  various  friction  models  on  the 
extrusion  process 


ALE 


•  Solid  mechanical  approaches  have  been  unable  to  deal  with 
the  enormous  deformations  (i.e.,  “flows”)  associated  with 
exjtrusion  problems 

•  The  reason  is  the  reliance  upon  the  Lagrangian  description 

•  Mesh  entanglement  requires  frequent  remeshing  throughout 
the  extrusion  process 


•  ALE  (arbitrary  Lagrangian-Eulerian  formulation)  has  been 
used  primarily  for  fluid  mechanical  phenomena  to  simulate 

free-surface  flows  and  fluid-structure  interaction 

•  One  of  the  fundamental  and  absolutely  essential 
technologies  for  extrusion  simulation 
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ALE  (cont.) 


•  Constitutive  integration 

•  Extremely  simple  for  typical  Newtonian  and  certain 
classes  of  non-Newtonian  constitutive  equations 

•  More  complicated  for  history-dependent  constitutive 
models,  nevertheless  widely  employed  in  chemical 
engineering  for  polymer  extrusion 

•  Stabilized  methods  widely  used  in  polymer  extrusion 


space-time  Formulations 

New,  potentially  superior  ALE  technology  has  been  devel¬ 
oped  in  recent  years  from  a  space- time  viewpoint 

Conservation  laws  automatically  enforced,  a  difficulty  with 
some  traditional  ALE  methods 

Comprehensive  mathematical  theory  exists 

A  rigorous  basis  for  a  posteriori  error  estimation  and  hence 
adaptive  refinement  strategies 

Efficiency  trade-offs  with  traditional  ALE 


Die  Flexibility 


•  It  is  essential  to  model  the  flexibility  and  heat  conduction 
behavior  of  the  container  and  die 

•  Requires  a  multiregion  thermomechanical  contact  algo¬ 
rithm  incorporating  a  flow-like  ALE  description  of  billet  and 
part  in  conjunction  with  a  classical  Lagrangian  description 
of  the  container  and  die 

•  Also  requires  a  multiphysics  solution  strategy 
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Region-based  Computational  Architecture 


•  The  system  is  composed  of  distinct  physical  regions 


•  Requires  a  computational  architecture  supporting  distinct 
physical  and  mathematical  models,  and  their  interaction 
through  interfaces 


•  Separate  regions  also  support  and  simplify  adaptivity  and 
instantiation  of  design  changes 


Parallel  Processing 


•  Rapid  progress  in  hardware  technology  and  practical 
necessity  suggest  extrusion  simulations  will  be  done  at  the 
desktop 

•  For  the  present,  these  calculations  will  be  CPU-intensive 
and  will  benefit  from  the  use  of  powerful  parallel  systems 

•  Two-level  parallel  architecture 

•  Top-level  adopts  the  message  passing  paradigm  on  a 
domain  decomposition  of  the  physical  problem 

•  Second  level  employs  the  data  parallel  paradigm 
within  the  subdomains 

•  Either  level  of  parallelism  can  be  absent 

•  Enables  efficient  performance  to  be  attained  on  a  broad 
array  of  parallel  architectures 

•  Operational  on  existing  platforms,  such  as  IBM  SP-2, 
CONVEX  Exemplar,  SGI  Power  Challenge,  Thinking 
Machines  CM-5,  etc. 


Equation  Solving  Technologies 


•  A  variety  of  parallelized  iterative  and  direct  equation 
solving  technologies  will  be  employed 


•  Continual  progress  is  being  made  in  this  area 


Visualization 


•  Modern,  three-dimensional  interactive  computer  graphics 
will  play  an  essential  role  in  the  visualization  of  process 
and  design  evolution 


•  Animation  and  archiving  capabilities  will  enable  a  physical 
understanding  to  be  attained  of  very  complex  phenomena 


Control  of  Process 


•  Code  architecture  needs  to  be  able  to  accommodate  user- 
friendly  programmable  process  control 


Leveraging  Existing  Technology 


Modeling  Environment 


AML  (Advanced  Modeling  Language)  from  Technosoft,  Inc. 

•  Feature-based  design  methodologies 

•  Parametric  modeling 

•  Knowledge-based  engineering  techniques 

•  Represent  designs  and  processes 

•  Manage  hierarchies  of  analysis  models 

•  Perform  model  reduction  and  feature  removal 


•  Effectively  utilize  design  information 


Mesh  Generation  and  Adaptivity 


Finite  Octree  from  RPI 


•  Interfaced  to  AML 


•  Automatic,  unstructured  mesh  generation 

•  Adaptive  remeshing 


Geometrical  Modeling  Kernel 


Shapes  from  XOX,  Inc. 

•  Interfaced  to  AML 

•  Provides  a  powerful  non-manifold  topology  environment  for 
building  analysis  models 
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Simulation  and  Visualization 


Spectrum  Multiphysics  Solver  and  Visualizer  from  Centric 
Engineering  Systems,  Inc. 

•  Extensive  multiphysics,  multiregion  parallel  simulation  and 
interactive  visualization  capabilities 

•  Integrated  with  AML,  Shapes  and  Finite  Octree 

•  Arbitrary  Lagrangian-Eulerian  meshing  capability  for  billet 
and  part 

•  Three-dimensional  free  surface  capability 

•  Extensible,  high-level  platform  for  research  and  develop¬ 


ment 


User  Interface  Toolkit  and  Optimization  Kernel 


ILOG  Views  from  ILOG 

•  Object-oriented,  graphical  user  interface  toolkit 

•  Customized  interfaces  for  end  users,  e.g.,  extruders,  die 
makers  and  designers 

Optimization  Kernel 

•  To  be  selected  when  the  research  reaches  a  sufficiently 
mature  stage 
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Summary 


•  Presently  the  design  of  the  die  for  extrusion  and  the 
selection  of  process  parameters  is  done  with  little  or  no 
automation 

•  This  process  relies  heavily  on  previous  experience  and  often 
involves  lengthy  and  costly  trial  and  error  cycles 

•  A  Simulation  Based  Design  system  for  extrusion  needs 
to  address  the  complete  thermomechanical  interaction  of 
the  billet,  die  and  part 

•  Successful  development  of  this  methodology  will  minimize 
the  process  design  cycle  and  lead  to  major  time  and  cost 
savings  in  product  development 


Automated  3-D  Mesh  Generation  for  Deformation 

Process  Design 

\  by  Prof.  Mark  S.  Shephard 
Scientific  Computation  Research  Center 
Rensselaer  Polytechnic  Institute,  Troy,  NY  12180--3590 

Abstract 

The  use  of  advanced  simulation  tools  by  metal  forming  process  designers 
is  hampered  by  the  complexity  of  developing  and  controlling  the 
numerical  analysis  models  required.  The  availability  of  automatic  mesh 
generation  capability  linked  directly  with  the  design  geometry  and 
reliable  analysis  methodologies  is  critical  to  the  elimination  of  this 
bottleneck.  This  presentation  focuses  on  techniques  being  developed  to 
automatically  generate  and  control  three-dimensional  finite  element 
meshes  for  use  in  deformation  process  design  driven  by  CAD  geometric 
models. 

An  automatic  mesh  generator  can  be  defined  as  an  algorithmic  procedure 
which  can  create,  under  program  control  and  without  user  input  or 
intervention,  a  valid  mesh,  or  geometric  triangulation,  for  geometric  models 
of  arbitrary  complexity.  The  development  of  algorithmic  procedures  for  the 
automatic  generation  of  finite  element  meshes  of  general  three-dimensional 
geometric  models  has  been  an  active  area  of  research  for  several  years  with 
procedures  which  can  work  with  CAD  representations  now  becoming 
available.  The  key  to  the  development  of  an  automatic  mesh  generator  is  an 
understanding  of  the  requirements  of  geometric  triangulation  and  to  ensure 
that  they  are  addressed  by  the  mesh  generator. 

The  'Finite  Octree'  procedure  has  been  developed  to  fully  meet  the 
requirements  of  an  automatic  mesh  generator.  Finite  Octree  generates  a 
mesh  in  a  two-step  discretization  process.  In  the  first  step  the  geometric 
domain  is  subdivided  into  a  set  of  discrete  cells  which  are  stored  in  a 
regular  tree  structure.  In  the  second  step  the  individual  cells  within  the 
Finite  Octree  are  discretized  into  finite  elements,  with  specific  care  to 
ensure  the  proper  matching  to  the  elements  in  the  neighboring  cells.  All 
of  the  individual  operations  performed  ensure  the  result  is  a  valid 
discretization  of  the  original  domain  at  that  point  in  the  process. 

An  important  aspect  of  an  automatic  mesh  generator  is  the  ability  to 
exercise  control  over  the  mesh  gradation  in  both  an  a  priori  and  a 
posteriori  sense.  Mesh  control  functions  which  are  spatially-based 
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provide  the  most  flexibility.  The  structure  of  the  Finite  Octree  mesh 
generator  provides  spatially-based  mesh  control  since  mesh  gradation  is 
controlled  by  the  tree  levels  -  with  each  octant  in  the  tree  corresponding  to 
a  known  portion  of  the  domain  geometry.  Automatic  mesh  generators  must 
interact  with  the  geometric  representation  of  the  domain  being  meshed. 

One  approach  to  provide  such  information  is  an  interface  that  directly 
employs  the  functionality  of  the  geometric  modeling  system  in  the  mesh 
generation  process.  Finite  Octree  employs  a  functional  interface  built  on 
the  geometric  modeling  system.  This  approach  is  consistent  with  the 
functional  interface  methods  of  the  CAM-I  Application  Interface 
specification. 

Finite  Octree's  functional  interface  uses  21  "point-wise"  geometric 
interrogations.  The  integration  of  Finite  Octree  with  a  new  geometric 
modeling  system  requires  the  creation  of  procedures  which  request  the 
geometric  modeler  to  perform  the  21  interrogations.  The  approach  has 
been  used  to  link  Finite  Octree  to  geometry  systems  such  as  Parasolids 
(which  is  embedded  in  the  Unigraphics  CAD  system),  Shapes  (a  product  of 
XOX  which  is  embedded  in  many  design  and  analysis  CAD  products)  and  ACIS 
(a  product  of  Spatial  Technologies  which  is  also  embedded  in  many  design 
and  analysis  CAD  products). 

The  availability  of  an  automatic  mesh  generator  does  not  ensure 
acceptable  accuracy  of  a  finite  element  analysis.  However,  the  careful 
combination  of  an  automatic  mesh  generator  coupled  with  an  adaptive 
analysis  system  can  accept  as  input  a  geometric  representation  of  the 
domain,  with  loads,  material  properties  and  boundary  conditions  applied  to 
the  geometric  model,  and  automatically  produce  solution  results  to  a 
prespecified  degree  of  accuracy. 

During  the  adaptive  analysis  of  an  evolving  geometry  problem,  factors 
that  may  indicate  a  desire  to  modify  the  mesh  include  the  need  to 
improve  the  distribution  of  elements  for  effective  control  of  the 
discretization  error  or  to  account  for  degradation  in  element  shapes  due 
to  large  deformations.  Tools  which  can  address  one  or  both  of  these 
needs  include  (i)  node  point  repositioning,  (ii)  local  refinement  and 
de-refinement  procedures,  and  (iii)  remeshing.  All  three  have  a  role  in 
an  efficient  metal  forming  simulation  procedure. 

Node  point  repositioning  procedures  change  the  mesh  by  repositioning 
the  grid  points  while  maintaining  the  given  mesh  topology.  Any  effective 


Adaptive  Lagrangian  Eulerian  (ALE)  procedure  must  have  a  node  point 
repositioning  procedure  for  the  positioning  of  the  grid  after  each  step. 
Procedures  for  positioning  of  the  grid  have  also  been  used  with  Lagrangian 
techniques  as  a  simple  remeshing  procedure.  Criteria  used  to  reposition 
nodes  can  consider  minimizing  the  discretization  error  (re-refinement), 
controlling  element  shapes,  or  a  combination  of  both. 

A  number  of  approaches  are  available  to  locally  refine  and  de-refine 
meshes  (i.e.,  hexagonal  3-D  mesh)  to  control  the  discretization.  When 
simplex  meshes  (i.e.,  tetrahedral  3-D  mesh)  are  used  there  are  a  number  of 
algorithmic  approaches  that  avoid  the  introduction  of  hanging  nodes 
(which  require  the  introduction  of  constraint  equations).  An  edge-based 
procedure,  with  local  swapping  provides  the  highest  degree  of  flexibility  in 
this  process.  Remeshing  procedures  have  been  used  in  two-dimensional 
metal  forming  simulations  to  both  control  element  shapes  and  to  reflect 
refinements  for  discretization  control.  Automatic  remeshing  has  the  same 
basic  requirements  as  the  automatic  generation  of  the  initial  mesh.  The 
added  functionalities  required  are  ensuring  the  automatic  mesh  generator 
can  reflect  the  spatially-based  mesh  gradations  defined  by  the  error 
estimators,  and  the  ability  to  update  the  geometric  representation  (as  needed 
for  meshing)  based  on  the  analysis  results  of  the  last  solution  step  before  a 
remesh. 
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AUTOMATED  3-D  MESH  GENERATION 
FOR  DEFORMATION  PROCESS  DESIGN 
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INTRODUCTION 

Application  of  advanced  simulation  tools  by  designers  is  hampered 
by  complexity  of  developing  and  controlling  the  required  meshes 
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AUTOMATIC  MESH  GENERATION 

Goal:  To  generate  a  valid  discretization  in  a  controlled 
manner  for  use  in  a  numerical  analysis  procedure 
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Essentially  same  effort  needed  for  either 
approach 


GEOMETRIC  TRIANGULATION 
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TOPOLOGICAL  COMPATIBILITY 
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GEOMETRIC  SIMILARITY 
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COMPATIBILITY  ASSURANCE 
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Can  use  general  algorithm  for  all  entitles,  however,  a 
tailored  approach  Introduces  computational  efficiency 


OCTREE-BASED  MESH  GENERATION 
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level  of  the  octree,  and  thus  the  cell  size 
through  the  object's  domain 
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DEPINITION  OF  FINITE  OCTREE 
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LOOP  BUILDING  ALGORITHM 
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-  Classification  and  type  of  loops  with  respect  to  model  face 

-  Mesh  and  model  adjacencies 

-  Parametric  spaces  of  model  faces  (when  available) 


Three  loops  of  edge-uses  and  one  vertex  loop  are  created 


OCTANT  TRIANGULATION 


ELIMINATION  OF  SMALL  SEGMEN 
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2.  Vertex-Edge  collapsing 

3.  Vertex-Face  collapsing 


Statistics  for  connecting  rod  modei 


boundary 

Two  step  procedure  developed  focusing  on  the 
boundary  and  poorly  shaped  elements 
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NODE  POINT  REPOSITIONING 


geometric  representation  and  shouid  be  able 
to  operate  on  non-manifold  geometric  domains 
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INTERACTIONS  BETWEEN  GEOMETRIC 
MODEL  AND  MESH  GENERATOR 


ent  platforms 

Time  required  to  write  the  interface  io  first  modeler  was 
months,  subsequent  interfaces  took  2  months 
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FINITE  OCTREE  INTERFACE 
TO  GEOMETRIC  MODELERS 


downstream  procedures;  no  mechanisms  formed 
no  loss  of  analysis  attributes  information 
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Elimination  of  Small  Mesh  Entities 


Size  of  feature 
1/5000**^  model  size 


1/5000’"  of  model  size  Worst  aspect  ratio  Worst  aspect  ratio 

=  261.916  =10.575 
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EXAMPLE  NON-MANIFOLD 
MODELS  AND  THEIR  MESHES 


BASIC  ISSUE  FOR  HEXAHEDRAL  ELEMENT 


FINITE  OCTREE  HEXAHEDRAL  SOLID  MESHES 


FINITE  QUADTREE/OCTREE 
AUTOMATED  FEM  SYSTEMS 
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with  refinement  graded  towards  boundaries 

•  Model  boundaries  are  locally  discretized 

•  Finite  elements  are  made  In  new  quadrants 
and  transition  zone 


Mesh  Generation  -  Local  Mesh  Updating 


RESIDUAL  BASED  ERROR  ESTIMATOR 
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and  for  quadratic  tetrahedra  the  shape  functions  used  are 

n:  =  ,  n;  -  ^^66^4 ,  n;  =  ,  n:  =  ^x6^|£4 


Adaptive  Results  -  Adaptively  Refined  Meshes 


Adaptive  Results  -  Accuracy 


BISECTION  BASED  MESH  PREDICTION 

Utilize  mesh  topoiogy  and  ciassification  knowiedge 
Exploit  hierarchy  of  mesh  entities:  examine 
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BISECTION  BASED  MESH  PREDICTION  (CONT.) 
Examine  to  further  define  distribution 
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Repeat  with  mT^mT®  (check  for  d( 
complete  definition 


IMPLEMENTATION  OF  MULTIPLE  LEVEL  REFINEMENT 

Loop  over  mesh  entities  in  spatial  order  d  =  0, 
Determine  max  new  size  requested  adjacent  to 


<  min/maxVA^  assign  additional  value  (/>[*  to 
Note:  >  coarsening 


MULTIPLE  LEVEL  IMPLEMENTATION  (CONT.) 
Specify  distribution  to  locai  remesh  as  refine  points 


MULTIPLE  LEVEL  IMPLEMENTATION  (CONT.) 


Necessary  model  entities  reinserted 
Elements  formed  in  areas  of  tree  modification 
Collapsing  and  smoothing  performed 


EXAMPLES 

Examine  convergence  (error  in  energy  norm) 

•  Asymptotic  rate  of  h-adaptive  convergence: 
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•  Exact  error  predicted  ^  ©  =  1 

SUN/SPARC  10,128Mbytes  real,1.7Gbytes  swap/temp 


3-D  ADAPTIVE  RESULTS  —  FLANGE  EXAMPLE 
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Efficiency  improved  if  singularity  degree  carefully 
approximated 


ALUMINIUM  EXTRUSION  DIE  SPECIFICATION 
Extrusion  of  Aluminium  is  done  using  flat-faced  dies 
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solid-shape  die  used  for  Aiuminium  extrusion 


MODEL  PARAMETERS  FOR  SOLID-SHAPE  DIES 
The  die  models  have  been  defined  using  the  AML  virtual 
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Complete  knowledge  of  the  interface  maintained 


MODEL  DEFINITION  WITH  AML 
VIRTUAL  GEOMETRIC  MODELER 

of  defining  die  models  has  been  automated 
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Output  from  AML  is  the  model  meshed  by  Finite  Octree 
Finite  Octree  queries  the  modeler  to  obtain  the  necessary 
information  about  the  model 


MESHING  DIE  MODELS  USING  FINITE  OCTREE 
Finite  Octree  mesh  generator  has  been  used  to  generate 
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has  been  used  to  generate  the  example  meshes 


ELEMENT  SIZE  SPECIFICATION  BASED  ON  MODEL  ENTITIES 
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Geometric  approximation  factor 


MESH  CONTROL  PARAMETERS 
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Well  suited  for  surfaces  with  variable 
e.g.,  free  from  surfaces 


MESH  CONTROL  PARAMETERS 
FOR  T  DIE  AND  MESH  STATISTICS 


Smallest  dihedral  angle _  4.594 

Largest  dihedral  angle  _  163.42 

Largest  edge/shortest  height  21 .892 


MODEL  AND  MESH  FOR  T  SECTION 


MESH  CONTROL  PARAMETERS 
FOR  L  DIE  AND  MESH  STATISTICS 


Smallest  dihedral  angle _ 9.969 

Largest  dihedral  angle _  148.455 

Largest  edge/shortest  height  10.760 


MODEL  AND  MESH  FOR  L  SECTION 


MESH  CONTROL  PARAMETERS  FOR 
CHANNEL  DIE  AND  MESH  STATISTICS 
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Smallest  dihedral  angle _  5.046 

Largest  dihedral  angle  _  153.164 

Largest  edge/shortest  height  19.436 


INPUT  TO  ANALYSIS  PROCEDURE 
Finite  Octree  generates  a  generic  neutrai  fiie  for  the 
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SOFTWARE  DEVELOPED  FOR 
WRITING  SPECTRUM  INPUT  FILE 
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TETRAHEDRAL  MESH  GENERATION  WITH 
MULTIPLE  ELEMENTS  THROUGH  THE  THICKNESS 

Rao  Garimella,  Mark  S.  Shephard 
Scientific  Computation  Research  Center 
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Applications: 

-  Multiple  elements  through  thin  sections 

-  Boundary  layer  generation 


Multiple  elements  through  thin  sections 

Define  “thin”  assuming  modei  discretization  available 
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model  boundary 


ILLUSTRATION  OF  THIN  SECTIONS 

Assume  3  elements  requested  through  thickness 


DETECTION  OF  LOCALLY  THIN  SECTIONS 

Forward  search  from  start  vertex  to  find  another  vertex 
classified  on  model  face 
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Average  discrete 
normal 


DETECTION  OF  THIN  SECTIONS  (CONTD.) 

Assume  elements  requested  through  thickness 

Forward  search: 
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Start  at  mesh  vertex  found  on  forward  search 

Go  to  adjacent  mesh  vertex  classified  on  model  face  that  is 
closest  to  start  vertex 


DETECTION  OF  THIN  SECTIONS  (CONTD.) 

Repeat  until  local  miminum  reached 
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Assuming  A//,  Is  the  number  of  boundary  mesh  vertices,  total 
computational  cost  of  finding  opposite  vertices  is  O(Nij) 


CREATING  MULTIPLE  ELEMENTS  THROUGH  THICKNESS 

-  Create  by  local  mesh  modification  (edge  spiit,  swap) 
procedures  made  as  efficient  possibie  for  this  application 
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SWAPPING  EDGES 

Swapping  of  edges  after  splitting  to 

-  align  edges  along  and  perpendicular  to  thickness 
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Transformation  from  “diagonal”  to  “zigzag”  triangulation 


SWAPPING  EDGES  (CONTD.) 

Swapping  sequence  attempts  to  convert  wedge  with  all 
“diagonal”  faces  to  all  “zigzag”  faces 
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EXAMPLE:  plate 


Number  of  elements  96  384  384 

Largest  Dihedral  Angle  150  118  96 


EXAMPLE:  Ring 


550  nodes 
2169  elements 


EXAMPLE:  ICCA  model  full  solid 
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EXAMPLE:  CENTRIC  DIE  MODEL 
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Future  work 

Capability  to  deal  with  multi-material  models 
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tetrahedronize  to  efficiently  get  final  mesh 
Speed  impro’/ements 


CLOSING  REMARKS 

Base  tools  to  support  automated  adaptive  analysis  in  deformation 
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3-D  geometry  procedures  when  remeshes  are  required 

Effective  transfer  of  solution  variables  when  remeshes  are 
required 


Alternative  Materials  &  Processes:  Incomplete 
Design  and  Multi-process  Optimization 

by  A.  Chemaly 
TechnoSoft,  Inc.,  and 
Visiting  Research  Scientist 
Wright  Laboratory 

Abstract 

Design  research  has  long  focused  on  geometric  aspects  of  coupling  product 
representation  and  analysis.  More  recently,  with  the  advent  of  feature-based 
design,  research  has  enabled  the  further  coupling  of  processing  knowledge. 

But  our  preoccupation  with  geometric  issues  has  deferred  the  response  to  this 
need  for  coupling  of  knowledge  about  shape,  material  and  process.  We  continue 
to  be  constrained  by  the  less  technical  issues  such  as  a  standard  set  of  features 
and  a  rigorous  explication  which  is  generic  and  unambiguous  across 
applications  and  vendors.  As  a  consequence,  progress  in  automating  process 
planning  and  material  selection  will  continue  to  be  beset  with  tedious  and  costly 
human  intervention. 

Materials  and  process  design  begins  with  part  specifications  and 
geometry,  and  requires  the  optimization  (in  terms  of  cost  and  time)  of  a 
sequence  that  describes  setups,  fixtures,  and  processing  details,  e.g.,  in  the 
machining  process  step:  operations,  tooling,  machine  data,  material 
properties,  together  with  step-by-step  process  visualization  . . .  beginning 
with  the  raw  material  and  ending  with  the  final  shape.  But  the  more 
challenging  goal  is  to  optimize  across  all  processing  steps,  while  comparing 
alternative  shapes,_materials,  and/or  processes.  While  the  spirit  of 
Intelligent  Knowledge-based  Engineering  (KBE)  is  to.facilitate  multi- 
disciplined  team  participation,  the  depth  and  complexity  of  the  design 
environment,  multiple  and  competing  materials  and  processes,  has  become 
an  obstacle  to  widespread  practice. 

The  current  trend  in  engineering  automation  systems  is  the  application 
of  artificial  intelligence  techniques  to  automate  the  production  of 
customized  machined  parts  by  integrating  the  product  design/constraints 
together  with  material  and  processing  knowledge.  The  intelligent  KBE  user 
interface  must  go  beyond  simple  part  geometry  and  tolerances  to  include 
functional  and  processing  constraints  along  with  material,  if  desired.  The 
part  geometry  is  typically  represented  in  terms  of  generic  features  such  as 
pocket,  slot,  boss,  etc.  along  with  user-defined  or  customized  features  called 
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profiles.  All  features  are  represented  in  terms  of  objects  that  can  also  have 
non-geometric  properties  to  aid  in  the  process  plan  generation.  For 
example,  the  raw  stock  object  would  have,  in  addition  to  its  geometric 
features,  property  attributes  such  as  hardness,  yield  strength,  impact 
strength,  operating  temperature. 

Each  feature  (geometric  and  non-geometric)  of  the  product  design  is 
defined  as  an  'intelligent'  object,  i.e.,  with  all  the  rules  required  for  its 
creation  and  interaction  with  other  objects.  The  system  enforces  design 
rules  such  as  over  or  under  constrained  designs  and  assists  in  the  process 
plan  generation  by  interrogating  the  processing  knowledge  base.  Unlike 
the  usual  KBE  (parametric  CAD/CAM)  system,  these  rules  need  not  only  be 
associated  with  the  part's  geometry,  they  could  also  relate  to  computational 
analysis  results  and/or  processing  resources.  The  part  geometry  is  just 
another  attribute  of  the  intelligent  KBE  system,  and  instead  of  simply 
evaluating  a  proposed  geometric  design,  an  intelligent  KBE  system  assists  in 
the  design  creation  process  through  three  primary  components:  1 ) 
objects/sub-objects,  2)  knowledge  bases,  and  3)  an  evolving  model. 

A  model  is  formed  from  objects  which  are  the  base  primitives  of  a  KBE 
system.  An  object  often  corresponds  to  a  real  world  entity.  For  example,  a 
crank-shaft  is  an  object  commonly  used  for  mechanical  design  of  engines. 
Part  objects  contain  ail  the  information  necessary  to  model  themselves. 
Objects  vary  in  intelligence.  Primitive  ones  merely  know  how  to  draw 
themselves.  More  intelligent  or  smart  objects,  use  the  knowledge  bases  to 
make  specific  design  decisions  regarding  size,  shape,  position, 
manufacturing  process,  and  so  forth.  Some  objects  even  decide  which  other 
objects  are  applicable  for  this  design.  Objects  can  also  be  formed  by 
grouping  together  other  objects  (sub-objects).  For  example,  an  engine 
object  may  be  a  grouping  of  a  crank-shaft,  a  piston  and  cylinder  assembly 
and  other  sub-objects  representing  various  lower  level  component  parts  of 
an  engine.  The  key  point  is  that  these  objec::  are  combined  together  to  form 
and  to  describe  a  part  model. 

Knowledge  bases  capture  a  set  of  object  associations  and  relations  to 
represent  the  domain  knowledge  of  a  product.  The  KBE  system  provides 
organized  access  to  this  knowledge  base  to  improve  the  design  process.  By 
automating  the  interaction  of  these  knowledge  bases,  the  design/evaluate/ 


revise  design  loops  can  be  minimized  or  eliminated.  The  interaction  is 
governed  by  the  relations  between  objects.  Each  object  includes,  and  thus  is 
constrained  by,  a  set  of  rules  which  are  applicable  to  the  object  which 
encapsulates  them.  The  KBE  engine  must  establish  the  order  of  dependency 
among  the  objects  which  therefore  organizes  the  interaction  between 
objects  and  their  respective  rules.  Note,  this  capability  of  the  KBE  engine  to 
establish  dependency  is  what  distinguishes  rule-based  systems  and  KBE 
systems. 

The  conceptual  part  model  describes  the  part  engineering  process  by 
combining  the  appropriate  part  objects.  It  describes  the  relationships 
between  major  and  minor  components  in  the  product.  The  conceptual  part 
model  is  often  represented  as  a  tree  structure  and  represents  the 
organizational  schematic  for  the  product  design  and  process  planning  cycle. 
The  basic  philosophy  behind  this  concept  is  that  a  design  is  composed  from 
features  as  basic  entities  instead  of  primitives.  In  addition  a  feature  based 
design  environment  supplies  easy  referencing  and  parametric  association 
for  feature  properties.  It  supports  cross-referencing  of  properties  within 
the  different  features  so  that  you  can  specify  relations  and  place  the  results 
into  other  properties. 

The  concurrent  engineering  philosophy  is  an  obvious  goal  related  to 
teamwork.  Alternative  design  specifications  are  required  when  a 
manufacturing  problem  is  encountered  in  a  product  that  cannot  be  built  or 
is  simply  too  expensive.  Obtaining  timely  and  quick  feedback  from  the 
different  engineering  disciplines  involved  in  the  design  process  is  the  key 
to  avoiding  rework  and  shortening  the  product  design  cycle.  Geometric 
reasoning  enables  the  system  to  analyze  part  features  to  automatically 
generate  parameters  for  process  planning  simulation  and  automation.  These 
parameters  enable  the  integration  of  the  part  design  activity  with 
manufacturing  and  inspection  planning. 

The  part  geometry,  design  specifications/properties,  and  the 
manufacturing/inspection  plans  are  maintained  and  stored  in  a  single 
object-oriented  part  model.  This  unique  object-oriented  part  model  enables 
engineers  from  various  areas  of  expertise  -  design,  analysis,  manufacturing, 
inspection  planning  -  to  interact  with  the  part  representation 
simultaneously.  As  a  consequence,  a  design  team  can  capitalize  on  the  near 
instantaneous  'what-if  benefits  of  concurrent  engineering. 


Implementation  of  a  design  process  that  efficiently  incorporates  the 
concepts  of  concurrent  engineering  is  the  real  challenge. 

KBE  systems  have  demand  driven  calculation  kernels.  Demand  driven 
calculation  based  systems  maintain  the  associativity  between  the  different 
objects  and  their  properties  relations.  Unlike  procedure  oriented  languages 
the  user  does  not  specify  the  order  of  computation,  objects  are  instantiated 
and  values  are  computed  only  when  they  are  demanded.  Demand-driven 
calculation  is  very  efficient  for  selective  analysis  of  large  models  to  inquire 
about  alternative  results  when  using  "what-if '  changes  to  ascertain 
different  property  values.  Properties  and  subobjects  are  not  instantiated 
until  they  are  demanded. 

In  object  oriented  modeling  systems  the  part  model  reflects  the  hierarchy 
of  the  design  intent.  The  part  model  structure  incorporates  the  relation 
among  all  the  various  parts  including  the  dependency  and  relation/ 
constraints  among  the  properties.  Dependencies  are  automatically  tracked 
and  are  used  to  trigger  the  constraints  to  enforce  the  relations.  Constraints 
and  dependencies  can  be  added  after  the  parts  and  subparts  are  created. 
Therein,  dependency  tracking  assists  in  augmenting  the  rules  and 
constraints  after  the  model  have  already  been  instantiated.  Dependency 
backtracking  automatically  verifies  rule  violations  and  notifies  properties  to 
recompute  according  to  constraints.  Modification  in  any  property 
automatically  updates  the  necessary  changes  in  the  entire  model  to  maintain 
consistency. 

Much  work  remains  in  feature-based  representation  architectures  to 
associate  geometric  design  with  process  and  material  knowledge.  This 
presentation  addresses  an  extension  to  features  which  is  predicated  on  an 
architecture  for  inteUigent  KBE  technology,  referred  to  as  INCOMPLETE  DESIGN, 
to  not  only  integrate  product,  process  and  materials,  but  to  also  organize  and 
apply  evolving  engineering  and  manufacturing  knowledge.  The  ultimate 
capability  is  to  capture  valued  expertise  and  avoid  relearning  costly  mistakes, 
while  capitalizing  on  past  successes  and  identifying  future  opportunities. 
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An  Evolving  Computational  Model  for  In-Situ 
Material  Aging  Assessment 

Yoh-Han  Pao 

Case  Western  Reserve  University, 
and 

Visiting  Research  Scientist 
Wright  Laboratory 

Abstract 

With  the  advent  of  self-improving  computational  models,  in-flight  self- 
assessment  of  a  composite  air  or  space  craft's  performance  limits  may 
become  a  reality.  Based  upon  the  air  or  space  craft's  specific  experiences 
( self-monitored  wear-out  history),  real-time  evolution  of  a  system 
performance  limits  model  will  continually  or  upon  request  inform  the  pilot 
or  maintenance  personnel  of  current  and  future  performance  limits  and 
structural  life  remainmg,  as  well  as  prediction  of  impending  failure.  The 
development  of  sensors  to  accomplish  this  task  has  been  an  ongoing  effort 
for  several  years  as  part  of  the  Air  Force  'Smart  Structures'  initiative 
Polymer,  metal  and  ceramic  composites  lend  themselves  well  to  this 
technology,  since  the  sensors  can  be  embedded  at  manufacture. 

Like  humans,  individual  aircraft  are  a  product  of  their  life  experiences  - 
a  unique  set  of  missions  characterized  by  pilot-specific  decisions  regarding 
speeds  and  maneuvers.  Also  akin  to  humans,  with  use  and  age,  the  strength 
of  materials  for  various  components  diminish  with  time.  Analogous  to  type 
A  humans,  military  aircraft  are  continually  driven  to  achieve  at  their  limits 
of  performance.  To  operate  safely  at  their  limits  these  systems  require  very 
accurate  assessment  of  the  strength  of  the  aircraft  and/or  individual 
components  or  to  be  informed  of  impending  failure  at  any  point  during  its 
useful  life. 

The  computational  method  must  have  the  potential  to,  in  real-time,  inform 
the  pilot  or  cockpit  crew  of  any  structural  limits,  either  due  to  progressive 
wear-out  or  instantaneous  damage  (projectile  impact,  etc)  and  the  magnitude 
of  those  limits.  This  is  useful  in  making  an  in-flight  advisory  of,  for 
example,  whether  the  aircraft  should  be  abandoned,  returned  immediately  to 
base,  complete  the  mission  with  restrictions,  or  complete  the  mission  with  no 
restrictions.  Once  the  mission  is  complete,  the  method  can  assess  whether 
the  structure  needs  immediate  repair,  can  be  flown  with  restrictions  (and 
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the  magnitude  of  any  restrictions)  to  depot  for  repair,  or  can  be  flown 
without  restrictions.  Beyond  these  safety  evaluations,  potential  cost  savings 
may  be  realized  with  this  'inspect  for  cause'  approach,  as  opposed  to  the 
'inspect  at  interval'  method  now  currently  used. 

In  principle,the  structural  integrity  and  functional  capability  of  structures 
made  of  composite  materials  can  be  monitored  with  the  use  of  distributed 
networks  of  embedded  or  imprinted  sensors.  The  challenge  is  to  have  the 
ability  to  make  sense  out  of  the  vast  arrays  of  information  items  afforded  by 
such  arrays.  This  presentation  addresses  neural  network  computing 
techniques  for  data  compression  and  'data  fusion'  which  facilitates  such  rapid 
meaningful  interpretation  of  bodies  of  in-situ  data.  Analysis  of  the  Donaldson 
and  Springer  approach  indicates  that  it  is  an  effective  one,  and  data  on  aging 
or  cumulative  damage  can  be  attained.  However  there  can  be  complex 
interaction  between  aging  mechanisms  on  both  microscopic  and  macroscopic 
scales.  This  paper  describes  network  computational  methods  for  'feature 
extraction'  and  for  'data  fusion'  which  can  be  applied  to  the  present  task  of 
building  an  evolving  model  of  the  state  of  a  composite  structure.  Illustrations 
of  such  application  to  benchmark  data  and  to  actual  sensor  data  are  described. 
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PRESENTATION  TOPICS 


1.  OBJECTIVE  OF  WORK 

2.  THE  MATERIALS  PERSPECTIVE 
Possible  to  Detect  Signs  of  Cumulative  Damage 
Interacting  Mechanisms 

3.  THE  STRUCTURES  PERSPECTIVE 
Interacting  Mechanisms 

Managing  and  Interpreting  Data  from  Net  of  Distributed  Sensors 

4.  THE  COMPUTATIONAL  MODEL 

The  Internal  Representation  Associative  Memory  and  Display 
Illustrative  Examples: 

Learning  Concepts  in  Reduced-Dimension  Space 
Interpreting  Sensor  Data 
Discovering  Associations 

5.  FUTURE  WORK 

Implement  Computational  Model  of  Distributed  Net  of  Sensors 
Validatewith  Real  Data 

Support  Development  and  Implementation  of  Methodology 
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Hgure  1 .  Illustrating  The  Idea  of  Monitoring  Load  Response  By  Coinparing  Ratio 
of  Actual  and  Reference  Loads,  At  Same  Value  of  Strain 


Figure  2  Illustrating  the  Difference  Between  Homogeneously  Stressed  Test 
Speckner>s  and  Non4lomogeneously  Stressed  Structures. 
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rigure  3  Illustrating  Ttie  Idea  of  Comparing  Observed  Response  to  that  of  Reference  and  Using 
Extrapolated  Ratio  to  Estimate  Ultimate  Failure  Load  for  ’Structure’. 

(Patterned  after  Figure  4  of  Donaldson  and  Springer  Reference 
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Figure  7  Predicting  Fracture  Failure  in  Random  Cyd'ic  Loading  and  In  Constant  Cydic  Loadfrig 


Applied  Load 


Figure  8  An  Exxagerated  Sketch  of  How  The  Predicted  Failure  Load  OecreMes  With  Every 
Cycle.  But  In  Practice  Extreme  Resolution  and  Accuracy  Would  be  Required. 
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Figure  1.  Exploring  The  Hyporthesis  That  The  Internal  Representation  Is  Meaningful 


Different  Internal  Representations  Obtained  with  Different  Initial  States 


A  Benchmark  Body  Of  Data 
Gasoline  Blending  Data 

"-“Bayesian  Control  in  Mixture  Models” 

By  L.  Mark  Berliner 

TECHNOMETRICS,  NOVEMBER  1987, Vol.  29,  No.4  P  455-460 
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Table  1:  Gasoline  Blending  Data 


xl  =  Butane 
x2  =  Isopenetane 
x3  =  Reformate 
x4  =  Cat  Cracked 
x5  =  Alkylate 

y  =  Research  Cetane  at  2.0  Grams  of  Lead/Gallon 


Reduced  Dimension  Representation  of  Gasoline  Blending  Data 
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Reduced  Dimension  Representation  of  Gasoline  Blending  Data 
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Learning  Classification  Rule 
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Exploring  The  Metric 
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No. 
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c 
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1 

AgGaS2 

0.28 

5.75722 

10.3036 
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2.55 

■a 

AgAlS2 

0.3 

5.73 

AgGaSe2 

0.27 

5.755 

10.28 

1.786 

1.8 

4 

CdSiAs2 

0.298 

5.884 

10.882 

1.849 

1.55 

5 

CdGeP2 

0.2839 

5.738 

10.765 

1.876 

1.72 

6 

AgAlTe2 

0.26 

6.296 

7 

CdGeAs2 

0.278 

5.9432 

11.2163 
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AgGaTe2 

0.26 
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1.1 
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1.916 
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CuAlSe2 

0.26 
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2.6 

12 
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0.25 
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1.959 
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6.09 

11.94 
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3.35 

■a 

ZnGeAs2 

0.25 
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0.75 

17 
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19 
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0.25 
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20 
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1.985 

1.24 

21 

CuTiSe2 

0.25 

5.832 
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1.994 

1.07 

22 

ZnSnAs2 

0.231 

5.851 

11.702 

iBi 

5.65 

11.3 

2.000 

1.66 

24 

ZnLnSe2 

11.614 

2.008 

0.95 

25 

CuLnS2 

0.2 

5.5228 

11.1321 

2.106 

1.54 

0.25 

5.555 

1.71 

Table  2:  Birefringence  Data 
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Pointers  to  Birefringence  Characteristics 
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Dimension 


No. 

t:l-5 

Fault 

1 

0.65190 

0.13019 

0.31398 

0.69901 

0.30067 

0.00000 

2 

0.27577 

0.56790 

0.24946 

0.61443 

0.70156 

1.00000 

3 

0.86528 

0.30303 

0.10538 

0.56716 

0.58797 

0.00000 

4 

0.15642 

0.83277 

0.58065 

0.37313 

0.58352 

1.00000 

5 

0.82369 

0.27834 

0.24731 

0.67413 

0.90200 

0.00000 

6 

0.35353 

0.67116 

0.16559 

0.65920 

0,82405 

1.00000 

7 

0.40958 

0.35241 

0.41290 

0.73881 

0.70601 

0.00000 

8 

0.35443 

0.33782 

0.55054 

0.70647 

0.71269 

1.00000 

9 

0.54702 

0.5735 

0.59355 

0.67413 

0.72606 

0.00000 

10 

0.34177 

0.60718 

0.79355 

0.79851 

0.64588 

1.00000 

0.83582 

0.74833 

0.00000 

12 

0.35353 

0.57800 

0.94409 

0.95025 

0.74610 

1.00000 

13 

0.47197 

0.32099 

0.36559 

0.58209 

0.52561 

0.00000 

14 

0.36528 

0.39843 

0.44731 

0.61940 

0.55457 

1.00000 

15 

0.44123 

0.29854 

0.34624 

0.57711 

0.55457 

0.00000 

16 

0.35805 

0.35354 

0.42150 

0.59701 

0.56793 

1.00000 

17 

0.49005 

0.32997 

0.41505 

0.72139 

0.67929 

0.00000 

18 

0.31284 

0.43547 

0.43656 

0.72388 

0.70601 

1.00000 

19 

0.43309 

0.31874 

0.39785 

0.71642 

0.73497 

0.00000 

20 

0.34991 

0.36251 

0.44946 

0.71144 

0.73051 

1.00000 

21 

0.46745 

0.26936 

0.40860 

0.69652 

0.72160 

0.00000 

0.37261 

0.42366 

0.70398 

0.70601 

1.00000 

fcl 

0.59042 

0.25253 

0.48602 

0.78358 

0.82628 

0.00000 

24 

0.38427 

0.37486 

0.48172 

0.79851 

0.80401 

1.00000 

25 

0.38156 

0.19753 

0.40645 

0.63930 

0.83296 

0.00000 

0.34810 

0.52189 

0.44516 

0.6890^ 

0.72160 

1.00000 

HI 

0.91134 

0.44301 

0.61194  1 

0.51225 

0.00000 

0.41863 

1.00000 

1.00000 

0.59453  1 

0.49220 

1.00000 

0.50723 

0.36364 

0.40645 

0.68159 

0.71715 

0.00000 

30 

0.34991 

0.47250 

0.45806 

0.70149 

0.70156 

1.00000 

31 

0.54069 

0.24691 

0.38279 

0.70647 

0.73051 

0.00000 

■?gl 

0.40404 

0.38710 

0.70149 

0.72383 

1.00000 

33 

0.41320 

0.32660 

0.41075 

0.68408 

0.71715 

0.00000 

34 

0.34991 

0.34007 

0.49247 

0.68906 

0.70379 

1.00000 

0.39873 

0.35354 

0.44516 

0.68906 

0.69710 

0.00000 

0.33906 

0.32323 

0.58065 

0.70149 

0.69710 

1.00000 

0.29747 

0.26824 

0.42366 

0.74378 

0.85746 

0.00000 

0.30561 

0.21886 

0.36129 

0.59950 

1.00000 

Table  3:  Typical  Sensor  Data 
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Reduced  Dimension  Representation  of  Sensor  Response 
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Materials  Design  Using  Pyramidal  Networks 

by  A.  Jackson 
‘  Wright  Laboratory 

Abstract 

Created  by  V.  Gladun  (Ukranian  Academy  of  Sciences)  in  the  late  sixties, 
pyramidal  nets  are  a  means  for  classifying  objects  into  a  directed  graph  that 
allows  rules  to  be  written  relating  the  objects  (instances)  and  their  attributes 
(data).  For  materials  design,  pyramidal  nets  allow  prediction  of  materials  and 
properties  to  be  done  using  empirical  data.  The  net  consists  of  attributes,  nodes 
representing  collections  of  attributes,  and  the  objects.  Once  the  net  is  created, 
special  nodes  can  be  identified  that  characterize  sets  of  objects  and  which 
allow  logical  conjunctions  and  disjunctions  to  be  formulated  that  provide 
insight  into  the  relationships  among  the  objects  or  the  attributes.  Two 
examples  are  presented.  The  first  deals  with  a  simple  set  of  properties  for 
materials  that  illustrate  the  generation  algorithm.  The  second  example  is  for  a 
set  of  chalcopyrite  compounds  and  illustrates  the  utility  of  the  pyramidal  net 
in  associating  properties  to  compounds  or  compounds  to  properties. 

Design  of  Materials  is  an  intensive  effort  requiring  contributions  from 
many  individuals  with  a  wide  range  of  backgrounds.  The  time  spent  designing 
a  new  material  can  be  considerable,  ranging  from  a  few  years  to  many.  In 
contrast  to  design  of  objects,  design  of  materials  has  not  received  as  much 
attention,  because  the  process  of  materials  design  is  rather  complex  and  the 
foundations  are  only  now  being  explored  and  articulated.  A  number  of 
approaches  are  being  considered,  including  neural  nets  of  various  varieties, 
genetic  algorithms  and  optimization  techniques,  and  various  forms  of  logic 
systems,  for  example  fuzzy  logic,  rough  sets,  and  pyramidal  nets,  the  subject  of 
this  presentation.  Created  by  V.  Gladun  in  the  late  sixties,  pyramidal  nets  are  a 
means  for  classifying  objects  into  a  directed  graph  that  allows  rules  to  be 
written  relating  the  objects  and  their  attributes.  For  materials  design, 
pyramidal  nets  allow  prediction  of  materials  and  properties  to  be  done  using 
empirical  data.  The  net  consists  of  attributes,  nodes  representing  collections  of 
attributes,  and  the  objects.  Once  the  net  is  created,  special  nodes  can  be 
identified  that  characterize  sets  of  objects  and  which  allow  logical 
conjunctions  and  disjunctions  to  be  formulated  that  provide  insight  into  the 
relationships  among  the  objects  or  the  attributes. 
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Data  for  10  Chalcopyrite  Compounds 

Compound  Density  c/aRatio  aParar  cPara  Gap 

AgAITe2  6.172  1.878  6.303  11.842  1.415 
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Implementation  of  the  Algorithm: 
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A  Rapid  Supervised  Learning  Neural  Network  for 
Function  Interpolation  and  Approximation 

by  C.'L.  Philip  Chen,  Senior  Member,  IEEE 
Department  of  Computer  Science  and  Engineering 
Wright  State  University,  Dayton,  OH  45435 
and  Visiting  Research  Scientist 
Wright  Laboratory,  WL/MLIM 
Wright-Patterson  Air  Force  Base,  OH  45433 

Abstract 

A  unique  implementation  of  the  Functional-Link  neural  network  is 
addressed  which  focuses  on  an  instant  learning  algorithm  to  rapidly 
compute  the  weights  for  supervised  learning.  For  an  n-dimensional,  N- 
pattern  training  set,  with  a  constant  bias,  a  maximum  of  N-r-1  enhancement 
nodes  are  required  to  learn  the  mapping  within  a  given  precision  (where  r 
is  the  rank,  usually  the  dimension,  of  the  input  training  set).  Based  on 
linear  algebra  theorems,  a  neural  network  which  is  similar  to  the 
architecture  of  the  Functional-Link  neural  network  is  proposed.  In  this  way 
the  learning  process  is  equivalent  to  the  process  of  finding  an  optimal 
weight  matrix.  Since  the  neural  network  is  a  flat  network,  the  optimal 
weights  can  be  solved  easily  through  linear  least  squares  computation.  For 
off-line  training,  the  proposed  network  and  algorithm  is  able  to  achieve 
"one-shot"  training  as  opposed  to  the  more  iterative  training  algorithms  in 
the  literature.  An  on-line  iterative  training  algorithm  is  also  presented. 
Similar  to  most  of  back  propagation  type  of  learning  algorithms,  the  given 
algorithm  also  interpolates  the  training  data.  To  eliminate  outlier  data 
which  may  appear  in  some  erroneous  training  data,  a  robust  weighted  least 
squares  method  is  proposed.  The  robust  weighted  least  squares  learning 
algorithm  can  eliminate  outlier  samples  and  the  algorithm  approximates  the 
training  data  rather  than  interpolates  them.  Several  experiments  show  very 
promising  results. 
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C.  L  Philip  Chen,  Senior  Member,  IEEE 

Department  of  Computer  Science  and  Engineering 
Wright  State  University,  Dayton,  OH  45435 
Visiting  Research  Scientist 
Wright  Laboratory,  WL/MLIM 
Wright-Patterson  Air  Force  Base,  OH  45433 

ABSTRACT 

This  talk  presents  a  unique  implementation  of  the  Functional-Link  neural  network  and  an 
instant  learning  algorithm  that  rapidly  decides  the  weights  of  the  supervised  learning.  For  an 
n  dimensional  V -pattern  training  set,  with  a  constant  bias,  a  maximum  of  N-r-l  enhance¬ 
ment  nodes  are  required  to  learn  the  mapping  within  a  given  precision  (where  r  is  the  rank, 
usually  the  dimension,  of  the  input  patterns).  Based  on  linear  algebra  theorems,  a  neural  net¬ 
work  which  is  similar  to  the  architecture  of  the  Functional-Link  neural  network  was  proposed. 
In  this  way  the  learning  process  is  equivalent  to  the  process  of  finding  an  optimal  weight 
matrix.  Since  the  neural  network  is  a  flat  network,  the  optimal  weight  can  be  solved  easily 
through  linear  least  squares  computation.  For  off-line  training,  the  proposed  network  and 
algorithm  is  able  to  achieve  "one-shot"  training  as  opposed  to  most  iterative  training  algo¬ 
rithms  in  the  literature.  An  on-line  iterative  training  algorithm  is  also  presented.  Similar  to 
most  of  back  propagation  type  of  learning  algorithms,  the  given  algorithm  also  interpolates  the 
training  data.  To  eliminate  outlier  data  which  may  appear  in  some  erroneous  training  data,  a 
robust  weighted  least  squares  method  is  proposed.  The  robust  weighted  least  squares  learning 
algorithm  can  eliminate  outlier  samples  and  the  algorithm  approximates  the  training  data 
rather  than  interpolates  them.  Several  experiments  show  very  promising  result. 
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Outlines: 


•  Least-Square  Neural  Networks 

□  A  variation  of  the  Functional-Link  Net 

•  One  Shot  Learning  for  off-line  training 

•  Iterative  Learning  for  on-line  training 

•  Function  Interpolation 

•  Function  Approximation 


•  Examples 
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Consider  a  mapping  from  an  input  space  ,  to  an  out¬ 
put  space, . 

Output  =  /  ( input ) 

•  Ann  dimensional  -pattern  input  training  set 

•  An  m  dimensional  A -pattern  output  training  set 

□  Statistical  Regression: 

Linear  and  Nonlinear  —  Closed  form  solution 

□  Neural  Networks  Approach 

Training  —  Analytic  solution 
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•  Neural  Networks  Approach 

Output  =  a  *  *  '  '  a  (a  (Input  *  Weights))  *  Weights) 
Weights 

where  a  is  a  nonlinear  function  such  as 


tanh  (x)  = 


sigmoid  (x)  = 


l-te 


-X 


,  or 


radial  basis  function  (x)  =  exp  (- 


llx-p.11^ 


2a^ 


) 


A  Neural  Net  Figure  here 
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Supervised  Learning  Approach 

•  Network  Topology: 

Multi-layer  vs  Functional-link 

•  Back  propagation  and  its  variants: 

Newton  (Hessian),  Conjugate  Gradient  etc. 

•  Disadvantage  of  BP: 

Slow  convergeness 
Retraining 
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Statistical  Regression  and  Least  Squares:  A  Review 

4  , 

ji  =  woi  +  w  jiXi  +  ■  •  •  +  [X„  +  ei 

y2  =  W02  +  ^12^1  +  •  ■  ■  +  W„2X„  +  £2  (2) 

“  ^Om  ^Im^l  ^nm^n 

where  y,  ,  x,-,  and  e,-  are  column  vectors, 

y;  =  I>i/>y2/-  ■  ■  ■  ^yNif^ 

X;  =  [Xu^2i>  ■  ■  ■  - 

£,'  —  [£l,-  >£2;  >  ]  • 

Expressed  in  the  matrix  form,  the  linear  model  of  input  and 
output  is: 


Y(Nxm)  -  A(Nx(n+l))W((n+l)xm)  +  ^(Nxm)  (3) 


In  practice,  the  least  squares  computation  is  to  minimize  the 
error 


<|)(iy,  £)  =  Z(x^wf-yj,)2  , 
k 


where  and  can  be  solved  by  making 
(|)(W,£)  a  minimum,  i.e.. 


3(|) 


=0  and 


dE 


=0. 


cP^/ 


In  general,  a  linear  family  of  basis  functions, 
"■  ,  gi;  such  that 


yj  = 


- 1  - 

(5) 


in  which  the  basis  functions  are  known  and  held  fixed. 
The  coefficients,  Cj ,  are  to  be  determined  according  to 
the  principle  of  least  squares. 


In  other  words,  we  want  to  minimize  the  expression 

N 

(t)(co,Ci,  •  •  •  ,c„)  =  2(Zc  g-(x,  )-y;)^ 

'  i 


The  necessary  conditions  for  the  minimum  is 


In  practice,  the  above  equations  may  be  difficult  to  solve: 

•  The  basis  functions  need  to  be  not  chosen  carefully 

•  The  basis  functions  should  be  linearly  independent 

•  Well  conditioned  for  numerical  computation 


<2 
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Nonlinear  regression  model: 

4  , 

y^-  =  f(x^-,9)  +  £/,  i  =  1,2,  '  '  '  ,n,  (8) 

where  0  is  a  vector  containing  p  parameters  and  p  <n 

•  Least  squares  formulation  for  finding  the  coefficients 

•  Much  more  complicated 


Linear  Algebra  Theorems 


Theorem  1: 


Consider  a  system  of  equations  AW  =  B  , 
and  denote  the  augmented  matrix  of  the  system 
A^  4  [  A  I  B  ].  Assume  that  W  is  an 
{n+\)xm  matrix. 


(1)  The  system  AW  =  B  has  a  solution  if  and 
only  if  rank  (A)  =  rank(Ag). 


(2)  The  system  AW  =  B  has  a  unique  solu¬ 
tion  if  and  only  if 

rank  (A)  =  n+1  =  rank(A^). 
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•  Supervised  Learning: 

Find  a  weight  matrix,  W,  given  an  A  and  a  F  matrix. 

□  Neural  Networks:  One  hidden-layer  architecture 

[  a{AW{)  ][  ^2  ]  =  B  4  o-Hy) 

where  is  the  weight  matrix  from  the  input  to  the 

hidden  nodes,  and  1^2  is  the  weight  matrix  from  the 
hidden  nodes  to  the  output  nodes. 

•  Increase  number  of  hidden  nodes,  hopefully, 

=>  increase  the  rank  of  the  matrix  [  a(Aiy  j)  ],  or 
rank(  [  o(AWi)  ]  )  =  rank(  [  oCAWi)  I  S  ]  ). 

=>  W  =A"‘  S  or  A+S 
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□  Neural  Networks:  Another  one  hidden-layer  architecture 
[  A  I  aiAW^)  ][  W2  ]  =  5  4  a-\y) 

Advantages: 

Rank  of  [  A  I  a(AWi)  ]  higher  than  [  aCAl^i)  ] 
More  independent  columns,  easy  to  find  a  solution. 
Trade-off: 

Additional  connections  from  input  to  output 


Input  layer 
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Functional-Link  Neural  Network 

•  Architecture 


target  values 
outputs 

(N+J+1)  parameters  for 
each  output  node 

enhanced  pattern 


•  Learning  Algorithm 


c25"7 
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Reiterate  the  Linear  Algebra  Theorem: 

Theorem  2:  Given  the  neural  network  architecture  mentioned 
above, 

(a)  the  network  needs  q-r  hidden  nodes  to  have 
a  mapping  form  A  to  5  if  and  only  if 

rank  (  [A  I  a(A  )  ]  )  =  rank  (  [  A  \  B  ])  =  q 

(b)  the  network  needs  maximum  N-r-\  hidden 
nodes  to  have  a  unique  mapping  form  A  to  5  if 
and  only  if 

rank  (  [A  I  a(A  ]  )  =  rank  (  [  A  \  B  ])  =  N. 


(3^  3 


-  14  - 

Adding  a  constant  bias  term 

•  Increase  the  rank  of  the  input  matrix  by  1 

•  Similar  to  regression  equation  (a  constant  term) 

•  Add  even  product  terms  to  function  approximation 

Consider  Taylor  expansion  of  the  nonlinear  function: 
Taylor  (o(A  )  )  =>  odd  terms 
•  If  add  a  constant  term 

Taylor  (Oy(w^^'“+  £  ,  where  is 

it=l 

the  constant  bias. 

•  Gives  us 


Uq  +  a 


bias 

jo 


n 


k=l 


n 

+  z 

k=\ 


X:uW 


bias  n3 


ik  '^jk 


r  + 


•  Further  expansion  on  the  above  equation  will  result 

even  product  terms  on  . 

•  The  same  reason  to  include  a  bias  in  the  a  function 

<y(AW'’‘‘^  +  P) 
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Algorithm  Rank-Expansion  with  Instant  Learning  (REIL). 

Input:  The  input  pattern  matrix,  A ,  the  output  matrix,  F,  pre¬ 
cision,  and  the  error  tolerance. 

Output:  The  weight  matrix,  W,  and  the  neural  network. 

Step  1.  Lop  off  output  elements,  7,  within  (-1,1)  according  to 
the  precision. 

Step  2.  Calculate  B  = 

Step  3.  Add  k  hidden  nodes  and  assign  random  weights, 
k  <N  -r  -1. 

Step  4.  Solve  weight,  W,  by  minimizing  I  lA  W  -  5  1 12. 

W  =  (A^A)-U^5  or  (W=A‘'B)  (9) 

Step  5.  If  mean-squared  error  is  not  satisfied,  add  additional 
nodes  and  go  to  Step  4;  otherwise,  stop. 

End  of  Algorithm  REIL 


c9<A:>0 


-  16  - 


Iterative  Learning 

•  Patterns  are  input 

^  T  ^  —  1 

•  Reduce  computing  time  for  (A  A ) 


•  Update  W(k+l)  based  on  W{k),  given  a  new 
new  b(k+l) 

•  Let  A  (^+1)  =  [A  (^) 

•  and  B(k+l)  be  partitioned  as  [B(k) 

a{k+l)  and  a 

Wik+l)  =  (AikfAik)  +  \A{kfB(k)  + 

(10) 

For  convenience  define: 

Pik)  =  [A(kfA{k)rK 

(11) 

Thus 

W(k)  =  P{k)A{k)B{k). 

(12) 

From  matrix  inversion  lemma,  we  know 

P(^+l)  =  [A{k)^ A{k) 

(13) 

=  P{k)  -  P{k)A(k)[l  +AikfP(k)A{k)]  'AikfPik). 


Q?UI 
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Using  Eqs.  (12)  and  (13)  in  Eq.  (10)  yields 

W(k+l)  =  W (k)  -  P (k)A  {k)[l  +  A{kf  P(k)A{k)]~^ 

[B{k+l)-A(k)Wik)]  (14) 

For  recursive  evaluation  of  P(k) 

•  P(k)  is  initialized  by  P(0)  =  (  1/ri)/, 

where  r|  is  a  positive  real  random  number  approaching 
to  zero  and  I  is  an  identity  matrix. 

•  Disadvantage: 

Very  difficult  to  find  P (k)  if  P(k)  =  [A{k)^ A(k)r^  is  singular 
A  better  method  is  being  discovered  now 


-  18  - 

Robustness  and  Function  Approximation  using 
Weighted  Least  Squares 

•  Problem  with  Weighted  Least  Squares: 

Interpolate  the  training  samples 

•  To  reduce  the  approximation  due  to  erroneous  collection  of 

training  samples 


•  To  eliminate  outlier  effect: 

An  outlier  datum  is  a  datum  whose  functional  value,  , 
is  somewhat  abnormal  (either  very  big  or  very  small) 
compared  with  rest  of  functional  values  of  other  points 
whose  abscissas  are  in  the  vicinity  of  . 

It  will  distort  the  fitted  value  and  prevent  the  approximation 
from  following  the  majority  of  the  training  samples. 


•  Approach: 


□  Define  residual  error:  =  y,  -  y^ 

where  y^-  is  the  desired  output  and 
y-  is  the  estimated  value. 
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□  Let  M  be  the  median  of  the  absolute  values 

of  the  residuals,  M  =  Median  Ir^  I 

□  The  robustness  weight  for  the  training  sample  point 

Xj  is  defined  to  be 


G(x,)  =  5(-^),  (18) 

where  p  is  coefficient  that  controls  the  value  of  the  bell-shape 
function  based  on  the  spreadness  of  the  residuals. 

The  median  absolute  residual,  M ,  is  a  measure  of  how  spread 
out  the  residuals  are. 

residual  is  small  (compared  with  pM)  robustness  weight  =>  1 
residual  is  large  =>  the  corresponding  weight  is  0 
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□  Let  M  be  the  median  of  the  absolute  values 

of  the  residuals,  M  =  Median  \  r^  I 

□  The  robustness  weight  for  the  training  sample  point 

Xj-  is  defined  to  be 


G(x^)  =  S{ 


(18) 


where  p  is  coefficient  that  controls  the  value  of  the  bell-shape 
function  based  on  the  spreadness  of  the  residuals. 

The  median  absolute  residual,  M,  is  a  measure  of  how  spread 
out  the  residuals  are. 

residual  is  small  (compared  with  pM)  robustness  weight  =>  1 
residual  is  large  =>  the  corresponding  weight  is  =>  0 
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The  weighted  least  squares  is  to  minimize  the  error 

<|)(lf ,  e)  =  IG  (x^  )(x^  w/  -  f,  (19) 

k 

Redefine 

A  ^  and  F  ^ 

Similarly, 

W  =  (a'^gA)-^A'^B.  (20) 

The  iterative  solution  for  weighted  least  squares: 

IF(A:+1)  =  (AikfGAik)  +  a[^iak+ir\A{kfGB{k)  + 

(21) 

and 

P{k)  =  [A{k)^GA{k)T^,  (22) 
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Simulation  Results 

•  Error  Measurement:  Mean-Squared-Error  (MSE) 

MSE  = 

where  \A  I  is  the  cardinality  of  A ,  is 

the  desired  output  and  y,-  is  the  estimated  value. 

Example  1:  Find  the  network  weight  for  the  given  problem  (  a 
constant  bias  "-1"  is  added  to  the  A  matrix): 
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Max-error : Fig.  4 


MSE:Fig.  5 


-24- 


Example  2:  Consider  a  nonlinear  classification  of  a  three  class 
problem. 


-  25  - 


Example  3:  Consider  a  mapping,  z  =  sin(x)sinO),  x,y  =0  to 

K. 

□  Training 

•  Forty-nine  data  each  with  0.5  interval 

•  10  hidden  nodes 

•  Average  training  10  seconds  in  an  Alpha  machine 

•  The  maximum  error  (between  the  desired  outputs  and 

the  actual  output)  ^  around  10“^. 

•  The  mean  squared  error  ^  around  10“^. 

□  Testing 

•  1089  data  points  each  with  0.1  equally  interval 

•  The  maximum  error  (between  the  desired  outputs  and 

the  actual  output)  ^  around  10“^. 

•  The  mean  squared  error  — >  around  10“^. 


£)1( 
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Example  4:  A  random  noise  signal  between  [  -ji,  |li  ]  is  added 
to  the  function  in  Example  3,  where  |i  =  10”^ 


Example  5:  A  combination  of  triangular  and  sine  waves. 
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Example  7:  Approximation  of  a  function  with  outliers 

•  The  weighted  least  squares  method  is  applied  ‘ 

•  p  =  6. 


Conclusion 


•  Rapid  Learning  aCnd  Iterative  Learning 

•  Least  Squares  and  Weighted  Least  Squares 

•  Linear  Algebra  Theorem 

•  Function  Interpolation  and  Approximation 

•  Simulation  Results 
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Constant-radius  Blending  of  Free-Form 
Parametric  Surfaces  and  New  Results  on 
Pythagorean-Hodograph  Curves 

Rida  Farouki 

University  of  Michigan,  Ann  Arbor,  MI 

Abstract 

A  key  requirement  in  the  geometrical  modeling  of  mechanical  components  is 
the  ability  to  smoothly  "blend"  two  surfaces.  This  problem  is  especially  difficult 
in  the  context  of  free-form  surfaces,  and  has  no  simple  solution  that  is 
compatible  with  the  standard  (Bezier/B-spline)  representation  schemes  of  CAD 
systems.  We  have  developed  a  numerical  method  to  construct  approximate 
constant-radius  blends  between  intersecting  parametric  surfaces  that  satisfy 
prescribed  tolerances.  The  procedure  is  based  on  using  a  "pseudo-arc-length" 
predictor-corrector  scheme  to  calculate  the  intersection  of  the  offsets  to  the 
surfaces,  and  the  corresponding  normal  projection  curves  on  those  surfaces. 

The  blend  is  parameterized  rationally  such  that  its  isoparametric  curves  are 
exact  circular  arcs  that  terminate  with  tangent  continuity  on  the  given 
surfaces,  and  are  centered  on  the  intersection  of  the  offsets.  Points  on  the 
"spine"  curve  of  the  blend  are  generated  uniformly  by  arc  length  along  it  (with 
adaptive  step  reduction  to  accommodate  for  regions  of  high  curvature)  by  the 
predictor-corrector  method.  This  circumvents  problems  encountered  with 
using  increments  in  one  of  the  surface  parameters  to  trace  the  curve. 
Interpolating  these  points,  and  the  points  of  tangency  to  the  given  surfaces,  by 
cubic  splines  gives  rise  to  a  tensor— product  blending  surface,  of  degree  2  in  one 
variable  and  9  in  the  other,  that  can  be  represented  exactly  in  rational  Bezier/B- 
spline  form. 

We  will  also  present  recent  results  concerning  Pythagorean-Hodograph 
(PH)  curves,  a  novel  class  of  free-form  parametric  curves  that  offer  important 
advantages  over  the  "ordinary"  polynomial/ rational  curves  used  in  CAD  systems 
(rational  offests,  exactiy-computable  arc  lengths,  better  curvature  distributions, 
etc.)  In  particular,  we  describe  methods  to  smoothly  interpolate  point  data  using 
PH  curves,  and  the  formulation  of  CNC  interpolators  for  driving  a  machine  tool 
or  other  device  along  a  curvilinear  path  according  to  a  feedrate  specified  in 
terms  of  the  intrinsic  geometry  (curvature,  arc  length,  etc.)  of  the  path. 


Footnotes: 

#  1  Joint  work  with  graduate  student  Ragnar  Sverrisson. 


Constant— radius  blends  for 
free— form  parametric  surfaces 


Rida  T.  Farouki  and  Ragnar  Sverrisson 

Department  of  Mechanical  Engineering  &  Applied  Mechanics, 
University  of  Michigan,  Ann  Arbor 


problem  statement 


given  intersecting  free-form  parametric 
surfaces  p{u,  v)  and  q(5,  t) 

construct  constant-radius  blend  surface 
S  that  meets  p(tt,  v)  and  q(5,  t)  with 
tangent  continuity 

S  is  (part  of)  a  “canal  surface”  —  i.e., 
the  envelope  of  a  sphere  of  fixed  radius  d 
whose  center  traverses  a  space  curve  C 

want  to  represent  the  blend  surface 
S  in  rational  Bezier/B-spline  format 
for  compatibility  with  CAD  systems 


geometrical  considerations 


•  the  spine  curve  or  “center-line”  C  for 
the  blend  is  the  intersection  of  the  offsets 

p>d{u,v)  =  p(^,'y)  +  dra{u,v) 
c\d{s,t)  =  q(5,t)  +  dji{s,t) 

at  distance  d  to  the  two  given  surfaces 

•  C  represents  the  trajectory  of  the  center 
of  a  ball  of  radius  d  that  rolls  in  contact 
with  p(u,  v)  and  q(5,  t) 

•  the  normal  projections^  V  and  Q,  of  C 
onto  p(it,  v)  and  q(5,  t)  are  “trim  curves” 
for  the  blend:  they  indicate  where  it  meets 
the  given  surfaces  with  tangent  continuity 


0^5/ 


mathematical  difficulties 


•  exact  closed-form  blends  available  only 
for  “simple”  surfaces  (e.g.,  natural  quadrics) 


•  want  blend  to  be  a  true  constant-radius 
canal  surface  of  relatively  low  degree 


•  the  canal  surface  corresponding  to  a 
general  spine  curve  C  is  not  rational 
(Pythagorean-hodograph  space  curves 
are  an  exception) 


•  want  the  trim  curves  V  and  Q  to  be 
isoparametric  loci  on  blend,  to  avoid  need 
for  “trimmed  surface”  representations 
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numerical  method 


“pseudo”  arc-length  predictor-corrector 
scheme  gives  (u,  v)  and  (5,  t)  values  that 
define  sequence  of  points  on  C  and  'P,  Q 


each  set  of  three  points  on  C,  P,  Q  defines 
a  circular-arc  section  of  the  blend  surface 


specially-formulated  rational  spline  surface 
interpolates  these  sections  such  that  all  its 
isoparametric  curves  are  circles  of  radius  d 
centered  on  C  and  terminating  on  P,  Q 


monitor  and  confine  approximation  errors 
beiow  prescribed  tolerance  by  an  adaptive 
step-size  procedure 


offset— surface  intersection 


offset  surfaces  at  distance  d 

Pd{u,v)  =  p{u,v)  +  dxri{u,v) 
(ld{s,t)  =  q{s,t)  +  dn{s,t) 


unit  normal  vectors  are 


m 


Pu  X  p„ 

)  u 

Pu  X  P„ 


X  qt 
q*  X  qt 


find  values  {uk,  Vk,  Sk,  tk}k=o  satisfying 


Rdiu,v,s,t)  =  Pdiu,v)  -  qd{s,t)  =  0 


points  on  C  : 
points  on  V  : 
points  on  Q  : 


rfc  =  Pdiuk,Vk)  =  qdisk,tk) 

Pk  Pdi^ki 

qfc  qd(SA:j  t/;) 


OFFSBT  SURFACE 


predictor  step 


let  C(A)  be  the  spine  curve 
parameterized  by  its  arc  length  A 

u(A),  v{X)  and  5(A),  t(A)  define 
spine  curve  in  parameter  space  of  surfaces 

C(A)  s  pd{u{X),v{\))  =  qrf(s(A),t(A)) 


differentiate  Rd(u,  v,s,t)  =  0  with  respect 


dv 


ds 


dt 


dR(i  du 

IT  ^  dA  dA  “  dA  ”  dA 


=  0 


i  •  Pd,t; 

j  * 

k  •  Pd,v 


-  i  •  qd,3 


-  j  •  <ld,s 

-  k  •  Cid,s 


• 

’  du  ‘ 

• 

- 1  •  <id,t 

dA 

1  •  Pd,u 

• 

ds 

• 

-  J  •  ^d,t 

dA 

J  *  P(i,u 

dt 

- 1 

1 

.  dA. 

- 1 

cu 

du 

dA 


3  equations  in  4  unknowns 


A  =  arc  length 


dC/dA  I  =  1 


dC 

dA 


du 


Av 


Vd,u  +  Pd,v 


du\^  dudv 

dXj  ■^^P''’“’P‘'’''dAdA 


1  quadratic  +  3  lin^^ar  equations  in  4  unknowns 
2  distinct  solutions  —  choose  one  that  gives 


dC 

dA 


dC 


dA 


it-i 


>  0 


hence  predictor  step  is 

du  du  ds  dt^ 
dA ’  dA ’  dA’  dA y  .  ^ 

motion  by  AA  along  tangent  to  C(A) 


{u.v,s.  :Li  -  {u,v,s,t)k  + 


correction  iterations 


first-order  Taylor  expansion  of  system 

Iid{u  +  V  +  Av,  s  +  As,  t  +  At)  =  0 

F{u  +  Au,  V  +  At>)  =  0 


“pseudo”  arc-length  condition 

dC 


F{u,v)  =  [C(A)-C(A,)] 


dA 


-  AA  =  0 


constrains  motion  to  be  perpendicular  to  tangent 


i  ■  Pd,u 

^  ‘  Pd,v 

-  i  •  q</,« 

-  i  qd,<  ' 

Au 

-i’Rd 

j  ■  Pd,u 

j  ’  Prf,t) 

-  j  • 

-j  •  ^d,t 

Av 

-yRd 

k  •  p<f,u 

k  •  Vd,v 

-  k  •  qj,. 

-kqd,< 

As 

-k.Rrf 

dC 

dC 

0 

0 

dA  ■ 

dA  ■  P'' ■” 

At 

-F 

solution  of  4  X  4  system  yields  correction  step 

=  («>  I'. i)Ui  +  iAu,Av,As,At) 
iterate  over  i  until  converged 
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blend  surface  construction 


INPUT  DATA 

points  ao, . . . ,  aAT  on  trim  curve  V 
points  bo, . . . ,  Uat  on  trim  curve  Q 
points  Co, . . . ,  cv  on  spine  curve  C 

each  data  set  (aj:,  bt,  Cjt)  defines  circular  arc 
of  radius  d  from  at,  to  bj,,  centered  on  c*, 

want  to  “skin”  sequence  of  circular  arcs 
with  rational  canal— surface  approximation 

rational  Bezier  form  of  conics: 

_  WaSi{l  -  uf  +  u;ss2(l  -  u)u  +  Wbhu^ 
~  Wa{l  -  uY  +  u;s2(l  -  u)u  +  WbV?' 

s  =  intersection  of  tangents  at  a  and  b 


a,  b,  c  — >  cubic  splines  a(i;),  b('?;),  c{v) 
interpolating  the  given  data 

“standard”  choice  for  weights 

Wa  =  Wh  =  I  and  Ws  =  cos^(j) 
does  not  give  rational  surface 

“non-standard”  weights 

Wa  =  Ws  =  cos^  10  and  Wb  =  I 
yields  rational  surface,  degree(ii,  v)  =  (2,9) 

only  an  “approximate”  canal  surface 
—  plane  of  circular  iosparametric  curves 
not  always  orthogonal  to  spine  curve  c{v) 
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Error  (degrees) 


further  work 


•  refine  adaptive  step-size  selection 
and  procedures  for  error  control 

•  reliable  method  to  select  starting 
points  for  predictor-corrector  scheme 


•  derive  expressions  for  Bezier/B-spline 
control  points  of  the  blend  surface 

•  find  lower-degree  representation 
for  canal  surface  approximations 


•  exact  canal  surface  formulations? 
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Parallel  curves  and 
Pythagorean  hodographs 


Rida  T.  Farouki 

IBM  Thomas  J.  Watson  Research  Center 
Yorktown  Heights,  NY  10598 
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Pythagorean-hodograph  (PH)  curves 


offset 


trimmed  offset  untrimmed  offset 


trimmed  offset  untrimmed  offset 


segments  of  untrimmed  offset  =  10,136 
segments  of  trimmed  offset  =  48 


2x105 


polynomial  hodograph 

x\t)  =  u{t) ,  y'{t)  =  v{t) 
maps  into  Pythagorean  hodograph 
x\t)  =  u\t)-v\t),  y\t)  =  2u{i)v{t) 
integration  gives  PH  curve  r(t) 

Pythagorean  2uv 
polynomial  triples 

vP  —  'iP 
{2uv,v‘^-u'^) 

offset  r(i)  +  d  i)  =  raaonal  cmye 

s{t)  =  I  \/x'^  +  dt  ^  J  +  v^  dt 
arc  length  =  polynomial  function  of  t 


30(j> 


real-time  CNC  interpolators 
for  Pythagorean-hodograph  curves 

specify  tool  feedrate  V  along  PH  curve 
as  function  of  arc  length  s  or  curvature  k 


CNC  “interpolator”  generates  reference  points  along  curve 
at  sampling  interval  At  according  to  specified  feedrate  V 

analytically  reducible  to  sequence  of  (monotone) 
polynomial  equations  in  the  following  cases  — 

V  inversely  proportional  to  k 
V  constant,  linear,  or  quadratic  in  s 
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feedrate  inversely  proportional  to  curvature 


feedrate 

#  cycles 

time 

constant 

210 

8.40  fis 

linear 
arc  length 

236 

9.44  fis 

quadratic 
arc  length 

261 

10.44  ^s 

inverse 

curvature 

178 

7.12  /is 

PH  quintics  are  more  "graceful"  than  ordinary  cubics 


?r  to 


complex  formulation 

real  and  imaginary  parts  of  the  square 
of  a  complex  polynomial  w(t)  =  u{t)  +  iv{t) 
form  a  Pythagorean  hodograph: 

r'(f)  =  x'{t)  +  i  y’{t)  =  w^{t) 

=  7i“(t)  —  v~{t)  +  i  2u{t)v{t) 


write  hodographs  of  adjacent  PH  quintics  as 
4{t)  =  [  at_i  (1  -  +  b*  2(1  -  t)t  +  at  ]  2 

4+iW  =  [  (1  -  t)^  +  bA:+i  2(1  -  t)t  +  a^+i  ] 

automatically  satisfy  condition 


2ajt  =  bjt  +  bjt+i 

continuity  at  nodes  A;  =  1, . . . ,  AT  —  1 

d"  3a|  +  3  (a^_i  +  a^)  b^  -I-  2bj^ 

=  15(p/t  -  Pit-i) 

PH  condition  for  arcs  k  —  I, . . .  ^  N 


3/q 


assume  end-derivatives  do  and  d^v  given: 
ao  =  do  and  =  d^v 


linearity  of  conditions 
allows  elimination  of  ai, . . . ,  a^v-i 


‘‘tridiagonal”  system 
of  N  quadratic  equations 
in  N  complex  unknowns 


=  17b?  +  3hl  +  12bib2 
+  14aobi  +  2aob2  +  12aQ 

-  60 (pi  -  po)  =0, 

fy(bi,...,biv)  =  3b?_i  +  27b?  +  3b?+i 

+  13bj_iby  +  13  bjbj.|.i  +  by_ibj-+i 

-  60  (py  -  Pj_i)  =0  for  ;  =  2, . . . ,  iV  -  1 , 

,  bjv)  =  17b^  +  3b^_i  +  12b;\^byv_i 
+  HaysTb^r  +  2aArbAr_i  +  12  a?^ 

-  60(pyv  -  P7V-l)  =  0 

defines  PH  quintic  spline 


Sts' 


homotopy  method 


system  to  be  solved: 

fi(bi, . . .  ,biv)  =  0,  . . . ,  fiv(bi, . . .  ,bAr)  =  0 

“initial”  system  with  trivial  solutions: 

gi(bi, . . . ,  biv)  =  0,  •  •  •  ?  ’  -I  ^n)  =  0 


continuously  deform”  initial  system  into  desired  system 
while  tracking  the  motion  of  all  solutions 

hy(bi, . . . ,  biv,  A)  =  A  fj(bi, . . . ,  bjv) 

+  (1  -  A)  z  gj(bi, . .  .,bAr)  =  0 

“homotopy  parameter”  A  G  [  0, 1  ] 


hi(bi, . . . ,  bAT,  A)  =  0,  hAr(bi, . . .  ,bAr,  A)  —  0 

homotopy  system  defines  solution  loci  in  x  R 
loci  are  non-singular  for  “almost  all”  z 


compare  with  “ordinary”  cubic  spline 
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A  Constraint  Driven  'Piping  Algorithm' 

Using  Cyclides 

Graduate  Students:  Seth  Allen  and  Vinod  Kumar 
‘  Faculty:  Deba  Dutta 
Mechanical  Engineering  Dept. 

Univ  of  Michigan,  Ann  Arbor,  MI 

Abstract 

Cyclides  are  (tubular)  algebraic  surfaces  of  degree  four.  They  are 
attractive  for  use  in  CAD  environments  for  several  reasons:  dual 
representations,  intuitive  defining  parameters,  etc.  Blending  applications  of 
cyclides  have  been  demonstrated  in  the  past.  In  this  presentation,  we  will 
describe  a  procedure  to  compose  cyclide  pieces  that  satisfy  designer  imposed 
spatial  constraints.  This  development  is  in  the  context  of  an  important 
application:  cable/wire  harness  design.  We  consider  the  wire  harness 
problem  and  propose  a  two-step  solution.  First,  generate  a  "configuration" 
that  staisfies  the  logic  diagram  specified  by  the  designer.  Next,  route  the 
cable  within  the  environment,  for  the  given  configuration.  The  constraint 
driven  cyclide  piecing  technique  is  used  in  the  second  step. 

Each  wire  bundle  is  represented  as  a  cyclide  tube  (which  in  turn  could  be 
composed  of  several  cyclide  pieces)  It  is  an  interactive  scheme  in  which  the 
harness  designer  "guides"  the  cables  (cyclide  pieces)  while  several 
constraints  are  automatically  satisfied.  Examples  of  constraints  include 
number  of  junction  points  in  the  configuration,  max/ min  bend  radii, 
max/ min  cable  dia,  max/ min  length  of  cable,  regions  in  the  environment  to 
be  avoided/included,  etc.  These  constraints  form  an  optimization-  problem  and 
using  Sequential  Quadratic  Programming,  we  obtain  the  collision  free 
geometry  of  each  cyclide  piece. 

The  first  step,  i.e.,  configuration  generation,  is  a  graph  theoretic  problem 
and  possible  approaches  include  Steniner  trees  and  algorithms  for 
(minimum)  spanning  trees.  Analogies  to  VLSI  and  network  synthesis  are 
being  explored  to  determine  a  viable  scheme  for  generating  a  good 
configuration.  Note,  for  a  given  logic  diagram,  there  exists  a  multitude  of 
configurations  and  currendy,  in  industry,  no  more  than  one  configuration  is 
considered.  In  terms  of  A¥  applications  and  relevance,  we  will  discuss  our 
problem  in  the  context  of  harness  design  problems  at  Lockheed  Martin  (air 
craft  design)  and  at  Phillips  Lab,  Kirtland  AFB  (wire  harness  issues  in  space 
satellite  design). 
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Work  related  to  Phillips  Lab,  Kirtland  AFB: 
Description: 
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There  are  wires  (approx.  82)  which  run  from  MWIR  to  the 
electronic  box  on  the  payload.  These  wires  bypass  the  VISS. 


Considerations  (contd.): 
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•  The  wires  adds  considerable  stiffness  to  the  VISS  which  needs  to 
be  mechanicallv  "soft". 
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The  Optimizer 

The  optimization  code  uses  a  standard  optimization  technique 
Sequential  Quadratic  Programming: 
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Optimization  Step 


a 

o 

•l-H 

4-> 

03 

N 


Cl. 

O 

C 

03 

O 

>4-> 

c 

QJ 

(Z1 

C/3 

•l-H 

CD 

c/3  j_, 

P 

QJ  rP 
~t->  QJ 


U 

CJ 

4-. 

C/3 

0) 


P 

CD 
C/3 

O 

o  ^ 

§  ^ 

>1-.  •<-> 

s  ^ 

^  .p 

O 

P  03 
03  t!! 


O) 

j:: 

H 


03 

C3JD 

O 

u 

Cl. 


P 

03 

V 

CO 


;h 

4-> 

CO 

p 

O 

u 


0) 

CO 

q; 

><-• 

CO 

0) 

Ph 

CO 

•  t^ 
•)_> 

03 

CO 


13 

0) 

4-4 

u 

CO 

(D 

X 

4-4 

u 

p 

CD 

P 

CO 

03 

P 

O 


CO 

O 

Cl. 

p— . 

03 

P 


a 

o 

4-^ 

P 

O  . 
•X3  <D 

CO  Xi 

r(  ^ 
P-i  (/) 

9 

4-4  CO 

CO  CC3 
^  O 

4-»  .P 
4-< 

CO  'rf 


CO 

O 

P. 


u 
CD 

P  'g 
O  P 

CJ  (D 


CD 

CJ 

CD 


U 

•  CD 

p 

O 

•1^ 

4^ 

03 

N 


P. 

O 

CD 


CO 

3 

03 

03 

> 

P 

CD 

CD 

(Xi 


6Qk 


the  vertex  of  the  tangent  cone  (3), 

the  center  of  the  specified  final  position  (3),  and 

the  radius  of  the  final  sphere  (1). 


e  Optimization  Problem: 
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Interactive  Step  (Contd.) 
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Constraints  Imposed  bv  the  Designer: 
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Stage  2:  Collision  Free  Cable  Routing 
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O  Start  @  Goal 


Approach  to  Routing  Problem 
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Ongoing  Work: 

•We  have  identified  that  the  configuration  problem  as  a  variant 
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Graph  theory  definitions  (contd.) 

•  The  spanning  tree  with  the  smallest  weight  in  a  weight  graph  is 
called  a  minimal  spanning  tree.  (MST) 

•  A  Steiner  Tree  (ST)  for  a  given  vertex  set  Z  (terminals  i.e. 


u 


O) 

CD 

0) 

rC 


t-t  CD 

a> 

2J  N  2' 


CD  C 

«  (3 

V  ^ 

0) 

C 


C 

o 

•43 

u 


g.) 

^  o;  s 

^  Sd  g 
0) 

T3 


"T3 

CD 


5:1. 

C 

CO 


^  «s  ^ 

i2  CD  ^ 

C  <2  cjo 

•  1-H  .C 

O  ’T3 
Q^S\  S 


c 

0) 


u 
0)  0) 


Jrt  ^ 

S  aj 

CD  > 


O 


aj  CD 

X  .tJ 

M-l  ^ 

•i-<i  GJ 

H  ^ 
CO 
W 


0) 

q; 
u 
C 
ca 

•M 

.  CD 

^  •1—1 


Ei; 


CO 

« 


<s 

u 


CO 


(« 

cu 

■13 

T3 

p 


ns  w 

a; 

a  ^ 

W  CD 


c 

(C 

T3 

<D 


H 

CO 


<D 
60 
T3 

CD  q; 

•I-H  ^ 

Q1  I'P 

g_t 

(D  O 

C 

O) 


X  CD 

60  CD 
Tj 
^  O 


^  « 
•S  ^ 

5i  ns 
>u 

CO  ^ 

cc  pg 

1  ^ 

33 

CD 

•  fH 

cu 
0) 
;-i 


CD 


H 

2 

CO 

w 


43  2 

.SP  6 

a;  "P 

^  I 

s  S 
s  * 

c  fr 

•rH  C/i 

6  w 

cn  o 

CD  ^ 
fC 

^  E— I  ^ 
^  ^  (S 

^  S 

c  ci  a 

•r*  15* 

<U  2 
a  ^ 

‘a  53 
55  *St 

*r*  ^ 

§  55 


CO 

0; 

43 


3dC^ 


Variants  of  the  Steiner  Tree  problem  —  Applications: 
•  Wire  routing  phase  in  physical  VLSI  design. 
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Graph  theory  definitions: 

•  Graph  G  =  (V,  E  )  consists  of  a  vertex  set  V  and  an  edge  set  E. 
Each  edge  is  identified  by  an  unordered  pair  of  two  vertices. 
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skeleton. 

•  If  G  is  weighted,  then  weight  of  a  spanning  tree  of  G  is  sum  of 
the  weights  of  all  branches  in  T. 


Overview  of  the  solution  strategy: 
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The  system  should  address  various  issues  like: 
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Case-specific  constraints  like  locations  of  clamps  etc. 


Proposed  Solution  Strategy: 

A  three  stage  methodology  is  followed: 


STAGE  3:  Do  a  feasibility  check  on  the  generated  solution  —  if  it 
satisfies  all  the  designer's  requirements. 


r  Philosophy: 
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A  better  approach  is  to  keep  the  user /designer  involved,  making 
it  an  interactive  system. 


Objective: 
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Generate  near  optimal  (i.e.  goo..^  solutions  rapidly  as  several 
iterations  might  be  required  before  a  final  solntiop  is  selected. 


Cable  Harness  Design 
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Unconstrained  cyclide  tube  composition  (contd.) 

•  The  cyclide  tube  is  represented  by  two  end  spheres  and  a 
tangent  direction  at  one  end.  The  tangent  at  the  other  end  is 
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constraints /requirements  and  avoid  collisions  with  parts. 


Solution  Strategy: 

STEP  1:  As  the  part  is  axisymmetric,  position  the  probe  (using  z  and 
(|))  so  that  its  mid-plane  lies  in  the  plane  containing  the  axis  of  the 
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constant.  Note  that  all  a  values  are  permissible  at  initial  position. 

STEP  5:  Divide  each  final  C-cell  (having  a  path  from  initial  C-cell)- 
finely  and  repeat  steps  3  &  4  to  exactly  reach  the  goal  point. 


Unconstrained  Cvclide  Tube  composition 


cyclide  tube 
composition 


EDDY  CURRENT  PROBE  SETUP 


Problem  Statement: 
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The  angle  a  is  set  as  a  degree  of  freedom  and  its  value  identifie 
by  the  algorithm. 


Degrees  of  Freedom: 
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Angle  between  sides 

Circle  on  cone/ cylinder  where  cyclides  meet  quadric 


Properties  of  Corner  Blends 
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Corner  must  be  symmetric;  i.e.,  angle  at  vertex  on  each  side  must 
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Constraint  driven  cyclide  tube  composition 
scheme  for  wire  harness  design 


